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Abstract 

What can we learn about quantum gravity from a simple toy model, without actually 
quantizing it? The toy model consists of a finite number of point particles, coupled to three 
dimensional Einstein gravity. It has finitely many physical degrees of freedom. These are 
basically the relative positions of the particles in spacetime and the conjugate momenta. The 
resulting reduced phase space is derived from Einstein gravity as a topological field theory. 
The crucial point is thereby that we do not make any a priori assumptions about this phase 
space, except that the dynamics of the gravitational field is defined by the Einstein Hilbert 
action. This already leads to some interesting features of the reduced phase space, such as a 
non-commutative structure of spacetime when the model is quantized. 



Outline and summary 



It is widely believed that in a quantum theory of gravity, the concept of a smooth spacetime 
manifold has to be replaced by something new, perhaps a kind of non- commutative spacetime, 
whatever this precisely means. In a non-commutative spacetime the position of a particle in space 
at a moment of time is not a well defined point on a manifold. There is some kind of uncertainty 
relation, which makes it impossible for the particle to be localized. The spacetime itself obtains a 
foamy structure at small length and time scales. It is no longer possible to measure, for example, 
the relative position of two objects with arbitrary precision. The length scale at which this is 
expected to happen is the Planck scale, which is about 4>i ~ lCT 33 cm. 

The purpose of the present article is to study a toy model, which is based on Einstein gravity 
in three spacetime dimensions. There are no local excitations of the gravitational field, and in 
particular no gravitational waves is a three dimensional spacetime [fil EH. But nevertheless it 
is possible to define non-trivial dynamical systems. The special feature of these toy models is 
that they have finitely many physical degrees of freedom, while at the level of Einstein gravity 
as a field theory they share all those feature with the full theory in higher dimensions which 
are interesting from the point of view of quantum gravity. These features are, in particular, 
the diffeomorphism group of the spacetime manifold, which appears as a gauge group, and the 
presence of a dimensionful coupling constant, which is Newton's constant G. 

There are two classes of such toy models. The topological models are those that do not include 
matter. Physical degrees of freedom only arise if the universe has a non-trivial topology [|3|, 
0fl. The particle models contain pointlike particles as matter sources. The physical degrees 
of freedom are then basically the relative positions and momenta of the particles in spacetime 
[§, There are many interesting and unexpected features of such particle models, some of them 
even related to time machines and black hole physics [[7], ||, |9|, In the context of quantum 
gravity, the basic idea behind the models is to use a pointlike particle to probe the structure of the 
spacetime at small length scales [jTT], [T2| [T3]] . If quantum gravity sets some principle limitation 
on the possible length scales that can be resolved, then this should be seen when the toy model 
is quantized. 

In three spacetime dimensions, Newton's constant G is an inverse mass or energy, in units 
where the velocity of light is equal to one. Its inverse defines the Planck mass Mpi, which is a 
classical constant in the sense that % is not involved in the definition. So, there is a natural energy 
scale in the classical theory of general relativity. Another special feature of three dimensional 
gravity is that it admits pointlike matter sources, which are not hidden in black holes. The 
gravitational field of a point particle with mass m is simply a cone with a deficit angle of SnGm 
[jl], ^ . The spacetime in the neighbourhood of the world line is the direct product of this cone 
with a real line. 

Due to this conical structure of the spacetime, there is a gravitational interaction between par- 
ticles, even though there are no local gravitational forces acting on the particles. Each world line 
is just a timelike or lightlike geodesies, thus a straight line in an otherwise flat spacetime. How- 
ever due to the non-trivial global behaviour of a geodesic on a cone, the particles are effectively 
attracted by each other when they pass each other. But actually, we are here not interested in 
these dynamical features the particles. There are various methods to solve the classical equations 
of motion, and to obtain an overview of all possible spacetimes |Q, ^ |14p. 
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What we are rather interested in are the kinematical features of the particle models. By this 
we mean, for example, the phase space structures, and in particular the way gravity influences 
the symplectic structure and the Poisson brackets. It is this what distinguishes gravity from other 
typical interactions that can be defined between point particles. It turns out that gravity not only 
enters the Hamiltonian, by adding some kind of interaction potential or whatsoever, but also 
the symplectic structure. When the model is quantized, this has an immediate influence on the 
commutation relations between certain operators, long before the actual dynamics is imposed 
by the Hamiltonian. So, we expect the principle effects of quantum gravity at show up at this 
kinematical level. 

It is therefore not of much help to know the classical solutions to the equations of motion, 
hence in this case Einstein's equations in a three dimensional spacetime. Instead, we have to 
define a proper action functional for the particle model, and derive from this the phase space 
structures. As we are dealing with a toy model that has no counterpart in the real world, we have 
to be very careful with this definition, if we want to learn something about quantum gravity and 
the way it influences the small scale structure of spacetime. In particular, we should not make 
any a priori assumptions, which might be motivated by the expected results. The only assumption 
that we are going to make is that the kinematical and dynamical features of the gravitational field 
are completely specified by the Einstein Hilbert action. 

The particle model is in fact well defined as a dynamical system at the level of Einstein gravity 
as a field theory on a fixed spacetime manifold. We shall show this in section |]. At this level, the 
definition of the particle model has more or less become standard [ |14| , [TJ, |I6| ]. However, there 
are then various different ways to proceed. The model is defined as a field theory, where the 
basic variables are the metric or the dreibein and the spin connection on a given spacetime. So 
it has an infinite dimensional configuration space, though only finitely many physical degrees of 
freedom. One can use the ADM formulation of general relativity [jlTD, to set up a Hamiltonian 
framework. But then one still has an infinite dimensional phase space. 

To this one can apply the usual phase space reduction methods. There are constraints and 
associated gauge symmetries, and these can be solved and divided out. What remains is a finite 
dimensional reduced phase space, which is basically spanned by the position and momentum co- 
ordinates of the particles. A possible way to derive this reduced phase space, which is motivated 
by the conventional methods for gauge field theories, is to impose a gauge condition directly on 
the basic fields, hence the metric. In the ADM formulation of general relativity, this means that 
the foliation of the spacetime is chosen in a particular way, introducing a globally defined abso- 
lute time and space. The gauge fixed metric on the space at a moment of time is then specified 
by finitely many variables, and roughly speaking they define the positions and velocities of the 
particles in space at that moment of time [14], [15], \Wj \. 

So far, this approach is completely straightforward. However, it turns out to be quite hard to 
return to physics after this mathematical framework has been set up. The reason is the following. 
By construction, the reduced phase space is covered by a global canonical coordinate chart. One 
has a Poisson bracket which is diagonal in the basic variables, and the Hamiltonian is some 
function of these variables. Formally, a state can be specified by fixing n points on a plane and 
their velocities, where n is the number of particles. However, this plane and the points on the 
plane have not very much to do with the real space in which the particles are moving. The real 
space at a moment of time is a conical surface with n tips, which is embedded into a locally flat 



2 



spacetime with n conical singularities. 

The actual geometry of this conical space is encoded in the locations and velocities of the 
points on the plane in a somewhat implicit way. Consider for example the following question. 
Given two particles, what is the distance of these particles in spacetime, at a fixed moment of 
time? The spacetime is locally flat, and thus a geodesic connecting two particles is just a straight 
line. There might be several such geodesies in a conical spacetime, but then we can choose one 
and ask for its length. Expressing this length as a phase space function is an interesting question, 
because it is the corresponding quantum operator which tells us something about the spacetime 
structure when the model is quantized. But unfortunately, in the gauge fixed formulation this 
simple quantity is a very complicated function of the basic phase space variables. 

To know the geodesic distance between two particles, we first have to know the spacetime 
metric in terms of the coordinates used, and to know this we have to take into account all particles 
and their relative motion. Hence, a very simple object, the flat spacetime metric, is encoded in 
a complicated way. The idea of this article is therefore to present an alternative phase space 
reduction, which leads to a different set of basic phase space variables. The price that we have 
to pay is that these coordinates do not longer provide a global chart on the reduced phase space. 
Instead, it will be covered by an atlas of finitely many local coordinate charts. But on the other 
hand, these coordinates will have an immediate geometric and thus physical interpretation. 

To explain this, it is convenient first to consider the limit G — » 0, where the gravitational 
interaction is switched off. In this case, the particle model reduces to a special-relativistic free 
particle model. The particles are moving freely in a flat three dimensional Minkowski space. 
The phase space of this model is spanned by the absolute positions and momenta of the particles 
with respect to the embedding Minkowski space, which serves as a reference frame. But it is also 
possible to introduce relative coordinates and momenta, and to define a phase space spanned by 
them. This will be shown in section []]. The idea is, in a sense, to use the free particle model as a 
much simpler and well known toy model for our real toy model, and to show the analogy to the 
ADM, or Hamiltonian framework of general relativity. 

We can think of the flat Minkowski space in which the particles are living as a spacetime 
manifold, which is foliated by a family of spacelike slices. They are labeled by a kind of ADM 
time coordinate, and the state at a moment of time is specified by the geometry of the respective 
slice. A typical such slice is shown in figure 0]. Its geometry is, on the other hand, specified by a 
set of relative position and momentum coordinates for the particles. They refer to the links of a 
triangulation. The time evolution of this geometry is provided by the mass shell constraints for 
the particles, in analogy of the Hamiltonian constraints of general relativity. 

So far, this is just a strange way to describe a very simply system on n free relativistic point 
particles. However, the remarkable feature of this description is that we only have to make a 
few modifications, or deformations, in order to describe the geometry of the spacetime of the 
interacting system. This will be the subject of section and the resulting phase space structures 
will be discussed in section |[ In a sense, we can switch on the gravitational interaction in a 
continuous way, using Newton's constant G as a deformation parameter. First we have to cut the 
polygons of the triangulation apart, as indicated in figure 0, then they are deformed, as shown in 
figure |3], and finally we can glue them together again. The result is a spacelike surface which is 
no longer globally embedded into Minkowski space. 

Instead, it is a conical surface with n tips. This is the space at a moment of time with the grav- 
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itational interaction between the particles switched on. The crucial point is thereby that the link 
variables, hence those phase space variables that define the geometry of the individual polygons, 
are still the same as before. In particular, they still represent the same physical quantities, for 
example the relative position of two particles in spacetime. This is very different to the gauge 
fixed method, where the phase space variables are also the usual position and momentum coor- 
dinates of the free particles in the limit G — > 0. But when gravity is switched on, the physical 
interpretation of the variables is lost. 

The only problem is that it is now no longer possible to use these geometric variables as 
global coordinates on the phase space. There are some principle obstacles, and therefore it 
is not possible to have at the same time a global coordinate chart on the phase space, and a 



straightforward geometric interpretation of the coordinates 018fl. So, we have to cover the phase 
space by an atlas of local coordinate charts. But as long as we are interested in the local features 
of the phase space, and the symplectic structure is a local object, we can stick to one of the charts 
and need not care about the global structure of the phase space. 

It is at this point, where the perhaps most interesting feature arise, namely the non-commutative 
structure of spacetime. Let us briefly explain what happens. For a pair of free particles in flat 
Minkowski space, we can introduce a vector z a (a = 0, 1, 2), representing the relative position 
of the particles in spacetime. Clearly, the Poisson brackets for the components of this vector are 
zero. Such a vector can still be defined when gravity is switched on, and it still has the same 
physical interpretation as a relative position of two particles in spacetime. However, now it turns 
out that the Poisson brackets and therefore the resulting quantum commutators are given by 

{z a , z b } = 8-kG e abc z c [z a , z b ] = 8m £ Pl e abc z c , 

where £pi = Gh is again the Planck length, and e abc is the Levi Civita tensor. It defines the 
structure constants for the three dimensional Lorentz group. So, we have in this case a well 
defined notion of what is meant by a non-commutative spacetime. The position coordinates of 
the particle no longer commute. And obviously this non-commutative structure disappears in the 
limit G -> 0. 

The same quantum commutator has also been found for a much simpler single particle model 
[|T3p. The only difference was that there the vector z a did not represent the relative position of 
two particles, but the absolute position of the particle with respect to a reference frame. Unfor- 
tunately, there was some ambiguity in the definition of the reference frame, and we had to make 
some additional assumptions, beyond the basic one about the Einstein Hilbert action. One could 
therefore argue a little bit about the result, and the suggestion was that these problems could be 
overcome when going over from a single to a multi particle model. 

And in fact, we shall here see that the non-commutative structure arises from the phase space 
reduction applied to the Einstein Hilbert action without any further assumption. The reason 
why the multi particle model is, in a sense, better defined than a single particle model has to do 
with the asymptotic structure of a three dimensional spacetime. To set up a proper Hamiltonian 
framework of general relativity, one has to impose a fall off condition on the metric at infinity. 
This is typically some kind of asymptotical flatness condition. Moreover, one has to introduce 
a reference frame at spatial infinity, which can be interpreted as the rest frame of some external 
observer. 

The phase space and the Hamiltonian are well defined only if such a reference frame is in- 
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eluded into the configuration space of the model [|T9|]. The Killing symmetries of the asymptotic 
metric are then of special interest. They represent the possible translations and rotations of the 
universe with respect to the reference frame. In the Hamiltonian framework, these are the rigid 
symmetries of the model, and the associated conserved charges represent the total momentum 
and angular momentum of the universe. Now, the special feature of Einstein gravity in three 
dimensions is the conical structure of spacetime, which also manifests itself at spatial infinity. 
The universe is not asymptotically flat, but asymptotically conical. 

The symmetry group of a cone, however, is smaller than the symmetry group of a flat 
Minkowski space. The only symmetries are time translations and spatial rotations, whereas 
boosts and spatial translations are not allowed. This restricted symmetry group is that of the cen- 
tre of mass frame of a system of point particles. It is therefore appropriate to choose the reference 
frame so that it coincides with the centre of mass frame. But on the other hand, this does not 
make sense for a single particle model, because there is then no physical degree of freedom left. 
And this was the reason for the various difficulties and ambiguities to arise, which were related 
to the definition of the reference frame. 

All these ambiguities disappear when we define the multi particle model so that the reference 
frame coincides with the centre of mass frame of the universe. In section ^ we shall show this 
explicitly, defining the model in the first order formalism of general relativity, using the dreibein 
and the spin connection. We'll see that the required boundary terms to be added to the Einstein 
Hilbert action, defining the reference frame at spatial infinity, are uniquely fixed without making 
any further assumptions, in contrast to the single particle model where this was not the case. But 
in all other aspect the multi particle model will be a straightforward generalization of the single 
particle model. 

In particular, the coupling of the particles to the gravitational field will be defined in the same 
way, by formally converting the matter degrees of freedom into topological degrees of freedom, 
and introducing generalized, group valued momenta of the particles. As all this is explained in 
detail in [Of]. So, we shall here only briefly sketch this. The actual subject of section |] is then to 
perform the phase space reduction, and to derive all those phase space structures which we have 
to define without further motivation in section J3L The logical order of the article is actually so 
that section |] comes first, then section |2| describes the model from the spacetime point of view, 
section ^| takes the phase space point of view, and finally in section [T] we take the limit G — > 0. 

There are then various ways to proceed, but this is not any more the subject of this article. 
It actually serves a rather technical purpose, namely to set up a proper classical Hamiltonian 
framework for the particle model, which is as general as possible, but also as simple as possible, 
in the sense that no further assumption enters the model beyond the basic one mentioned above 
about the Einstein Hilbert action. The results can therefore serve as a starting point for a further 
analysis after some simplifications are made. For example, the simplest case is the Kepler system 
with n = 2 particles. This can be solve completely both at the classical level, deriving the 
trajectories in the phase space, and at the quantum level, solving the Schrodinger, or actually a 
generalized Klein Gordon equation [ |20| , ^TJ] . 

And for this particular example, we find indeed some interesting features of quantum gravity. 
It is no longer possible to localize the particles in space, or to bring the particles closer to each 
other than a certain distance, which is of the order of the Planck length. All this can be seen 
very explicitly and clearly when the model is quantized canonically. It is even possible to derive 
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the wave functions explicitly and see how the energy eigenstates of look like, when expressed in 
terms of physical variables. Hence the states can directly be interpreted as probability amplitudes, 
and we see explicitly how quantum gravity restrict the ability to probe small length scales. 

Another way to proceed further might be to consider the low coupling limit, in which case 
it could be possible to solve and quantize the n particle system exactly, or to consider a closed 
universe with a few particles as a new kind of cosmological toy model. But whatever one does in 
this direction, solving always means more then just writing down the solutions to the equations 
of motion or the quantum states in terms of some variables. This tells us almost nothing about 
the physics. To learn something about quantum gravity one has to give the solutions, hence the 
classical trajectories or the quantum states, in terms of physical variables. And for this purpose 
the Hamiltonian framework presented here might be useful. 



1 Free particles 

Consider a very simple system of n > 2 free relativistic point particles, living in a three dimen- 
sional, flat Minkowski space. As a vector space, we identify this with the spinor representation 
sl(2) of the three dimensional Lorentz algebra. The position variables x n = x^ a and the mo- 
mentum vectors p w = p^ a of the particles are represented as traceless 2x2 matrices, with 7 a 
(a = 0, 1, 2) being the orthonormal basis of the usual gamma matrices ( |A.1| ). A collection of 



definitions and formulas regarding the vector and matrix notation can be found in the appendix. 
The index n is used to label the individual particles. 

We have a 6n dimensional phase space which is spanned by the variables x n and p^. The 
symplectic potential and the resulting Poisson brackets are 

6 = ]T |Tr(p^d^) {p«,4} = rf h . (1.1) 

7T 

Note that we sometimes switch between the matrix and vector notation, and use whatever is 
more convenient. The Hamiltonian is a linear combination of the mass shell constraints C^, with 
Lagrange multipliers £ n E M. as coefficients, 

H = Y,^C*, C^ = \Ti(p^) + \m^Q. (1.2) 

7T 

The physical phase space is the subspace defined by the mass shell constraints and the positive 
energy conditions, 

\Tr{p%) = -ml pi = |Tr(p w7 °) > 0. (1.3) 

The mass shell constraints are first class constraints, and the associated gauge symmetries are 
the reparameterizations of the world lines as functions of a common, unphysical time parameter 
t. The gauge freedom corresponds to the freedom to choose the multipliers in the time evolution 
equations, 

Pn = {H,Pn} = 0, X n = {H, X^} = CttPtt- (1.4) 

The dot denotes the derivative with respect to t. The free particle system is a simple toy model 
for the ADM formulation of general relativity. The coordinate t on the world lines is the glob- 
ally defined, but unphysical ADM time. The mass shell constraints C n are the Hamiltonian 
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constraints. The multipliers ( w represent the lapse function. And the gauge symmetries, the repa- 
rameterizations of the world lines, are the spacetime diffeomorphisms. To make this analogy 
even closer, we impose the following two restrictions on the phase space variables. The first is 
a gauge restriction. It is analogous to the restriction to spacelike foliations in general relativity. 
The second restriction is a kind of asymptotical flatness condition, although at the moment this 
is not immediately obvious. 

The gauge restriction is as follows. At each moment of ADM time t, the particles must be 
located on a common spacelike surface in the embedding Minkowski space. In other words, the 
parameterization of the world lines is induced by a foliation of Minkowski space by a family of 
spacelike slices. For simplicity, let us further assume that the states where the particles are at the 
same point in space are excluded from the phase space. As long as we do not study the dynamics, 
that is the classical trajectories, this is not going to be a problem. The gauge condition is then 
requires all relative position vectors to be spacelike. Hence, 

spacelike for all n 1 ,n 2 . (1.5) 

With this restriction, the phase space is no longer a vector space. But it is still a 6n dimensional 
manifold, because the set of all spacelike vectors is an open subset of Minkowski space. The 
number of independent gauge symmetries in not affected by this restriction. 

The second restriction has to do with the rigid symmetries. The positions x n and momenta 
p v are defined with respect to a reference frame. We can think of it as the rest frame of some 
external observer. The rigid symmetries are the translations and Lorentz rotations of the particles 
with respect to this reference frame. The associated Noether charges are the components of the 
total momentum and angular momentum vector, 

P = Ep^> J = IS \Pn,x n ]. (1.6) 

7T 7T 

The reference frame coincides with the centre of mass frame if 

P = M 7o , J = S lQ , M,5 6l, M>0. (1.7) 

This is the second restriction that we impose on the phase space variables. In the centre of mass 
frame, there exists an absolute time and an absolute space. The symmetry group is restricted to 
the two dimensional subgroup of time translations and spatial rotations. A plane orthogonal to 
the 70-axis represents an instant of absolute time. The 7 -axis is the world line of & fictitious 
centre of mass particle. The total energy M is the mass, and the total angular momentum S is the 
spin of this particle. These are the conserved charges associated with the remaining symmetries. 

Except for one very special situation, a centre of mass frame always exists. If all particles 
are massless, and if they move with the velocity of light into the same spatial direction, then the 
total momentum is lightlike. These states are excluded. The dimension of the restricted phase 
space 6n — 4, since only the four spatial components of the equations Q1.7| ) are restrictions on 
the phase space variables. The time components are the definitions of the charges M and S. The 
gauge symmetries are not affected. We still have the freedom to reparameterize the world lines 
independently, within the range that is allowed by fll.5|). 

From the phase space point of view, the restriction to the centre of mass frame provides two 
pairs of second class constraints, whereas the mass shell constraints are still first class. For a 
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Figure 1 : The triangulated ADM surface of the free particle system. The particles are located 
at the positions x n in the embedding Minkowski space, which defines the centre of mass 
frame and also the reference frame. The momentum vectors p n are the tangent vectors of the 
world lines attached to these points. The surface is cut along the links A of a graph r, and 
divided into a collection 77 of simply connected polygons A. 



proper Hamiltonian description of the restricted system, we either have to replace the Poisson 
brackets with Dirac brackets, or we have to eliminate the second class constraints, by going over 
to a new set of independent phase space variables. This is what we are going to do now. We 
are looking for a convenient way to parameterize the solutions to the equations P = M70 and 
J = S~fo, which can later be generalized to describe the interacting particles. 

Triangulations 

Consider a foliation of the embedding Minkowski space by a family of spacelike slices, labeled 
by an ADM time coordinate t, and so that the parameterization of the world lines is induced by 
this foliation. A typical slice is shown in figure [T]. Let us call it the ADM surface at a moment 
of time t. Note that we have to keep the ADM time t apart from the absolute time defined by 
the 70-coordinate in the centre of mass frame. The ADM surface need not be an spatial plane 
in Minkowski space. For the time being, it can be any spacelike surface. The positions x n of 
the particles at the given moment of time are the intersections of the world lines with the ADM 
surface, and the momenta p w are the tangent vectors of the world lines attached to these points. 

A triangulation is defined by a collection r of oriented links. Every link X E T represents 
a geodesic in Minkowski space. It either connects two particles or extends from a particle to 
infinity. The particles are thus the vertices of a graph, and there is an additional, special vertex 
at infinity. The graph must be connected, hence there must be at least one link attached to 
every vertex, and it must be possible to go from every particle to every other along the links. 
Furthermore, the graph must be spacelike, which means that the ADM surface can be smoothly 
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deformed, so that all the links are contained in the surface. And finally, let us also require that 
all the links extending to infinity are spatial half lines. A spatial half line is a geodesic which is 
orthogonal to the 7 -axis, and extends from some point to infinity. 

If this is the case, then the graph r divides the ADM surface into a collection IJ of simply 
connected polygons. To define the various relations between the vertices, the links, and the 
polygons, we introduce the following notation. With each link A G r, we also have the reversed 
link —A G r. It represents the same geodesic with opposite orientation. We distinguish between 
internal and external links. An internal link A begins at a particle 7r_ A , and it ends at a particle 
ir\. For the reversed link — A, these two particles are interchanged, which explains the notation. 
The subset of all internal links is called T C r. 

It is sometimes useful to assign a preferred orientation to each internal link. For this purpose, 
we split the set r = T + U T_ into two disjoint subsets, so that A G i~+ implies —A G -T_ and 
vice versa. Since a priori there is no such preferred orientation, this decomposition is arbitrary. 
Whenever we use it we have to make sure that the results are independent of the particular 
decomposition. The external links already have a preferred orientation. We define C T to be 
the set of all external links oriented towards infinity, and r_oo is the set of all reversed external 
links oriented towards the particles. Again, it is so that rj E implies — rj G fL^ and vice 
versa. 

An external link rj G begins at some particle 7T_ V , but it does not end at any particle, and 
vice versa for the reversed external link — rj G We use the symbol A to denote a link 

in general or an internal link, whereas the symbol 77 always denotes an external link. A useful 
identity is the following decomposition of the graph r into four disjoint subsets, 

r = r + ur_ur 0O ur„ oo , r = r + ur_. (1.8) 

If 71 is some vertex, then r n = {A | n\ = n} is the set of all links ending at it, and r_ n = 
{A 1 7r_ A = 7r} is the set of all links beginning at ix. For the special vertex at infinity, this agrees 
with the definition of and fl^ above. Since every link begins and ends at some vertex, 
including the one at infinity, we have the following alternative disjoint decomposition of the 
graph, 

r = r ro u |J r n = u (J r_ n . (1.9) 

It is useful to assign a cyclic ordering to each of these sets. Let n be a particle, and consider a 
loop surrounding this particle in clockwise direction. The ordering in which the links cross this 
loop defines the cyclic ordering of the set r w . If A G is a link ending at it, then the next 
following link in clockwise direction is called the successor of A, which is denoted by A' G r^. 

If there is only one link attached of a vertex, then we have A' = A, otherwise this obviously 
defines a cyclic ordering. The same applies to the set of external links, but for consistency 
this is ordered in counter clockwise direction. Thus every external link i) G 4 has a unique 
successor rj G F^, which is the next external link that we cross when we walk around a circle at 
infinity in counter clockwise direction. Since every link is contained in exactly one of the subsets 
r n or Too, every link also has a unique successor, which is an element of the same subset. 

Each pair of links A and —A represents the common boundary of two polygons A A and A_ A . 
We define A A to be the polygon that lies to the left of A, and consequently A_ A lies to the 
right of A. For each polygon A G 77, the set 7a = {A | A A = A} represents the boundary of 
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the polygon, traversed in counter clockwise direction. A link A G /a is called an edge of the 
polygon A. Since every oriented link is an edge of exactly one polygon, we have yet another 
disjoint decomposition of the graph, namely 

r= U r A . (l.io) 

Aei7 

There is relation between the cyclically ordered sets r n and J^, and the sets i~A of edges. If 
A, A' G r w are two successive links ending at a vertex tt, which might also be the one at infinity, 
then A, —A' G i~A are two successive edges of some polygon A. For A, A' G r n we always have 
A A = A_ A ,. 

We distinguish between compact polygons A G ilo, and non-compact polygons A G /Too, so 
that II = n U iloo. A compact polygon A G 77o is bounded by at least three internal edges 
A G /a C i~o. A non-compact polygon has exactly one external edge r] G i~A H Too oriented 
towards infinity, exactly one reversed external edge — rf G Ja. H /"-oo, and one or more internal 
edges A G r A H .T . And the links r], rj' G i"^ are thereby always two successive external links. 
Consequently, there are as many non-compact polygons as there are external links. 

There is also a relation between the number of compact polygons and the number of internal 
links. It follows from Euler's polyhedron formula, and it involves the number of particles. Let 
£o = = |r_| be the number of internal links, so that each link is only counted with one 
orientation, and let = \Tool = |r_oo| be the number of external links. Similarly, p = |77 | 
is the number of compact polygons, and p^ = |77oo| is the number of non-compact polygons. 
Then we have the following relations, 

£ -pQ = n-l, 4o = poo. (1.11) 

They basically define the topology of the surface which is triangulated. We may release these 
conditions later on when we consider the coupled system, considering more general spatial 
topologies. But for the free particles the ADM surface always has the topology of R 2 . 

Relative momenta 

Let us now come back to the phase space of the particles. The primary problem was to solve the 
second class constraints (|1.7|). Let us first consider the equation P = Mj , which only involves 
the momenta of the particles. To find the general solution, we introduce a relative momentum 
vector q\ G sl(2) assigned to every link A G T, so that 

Ptt = ^<?a, q\ = ~q\- (1-12) 

The momentum p n of a particle n is the sum of the relative momenta q\ of all links ending at n. 
It is useful to think of the vectors <7a as currents flowing through the links. The second equation 
then states that the current flowing through the reversed link is the reversed current. And the first 
equation states that the momentum of a particle it is the total current that flows out of the graph 
at the vertex it. Using this picture, and the fact that the graph r is connected, it is easy to show 
that, for a given set of momenta p n , it is always possible to find a suitable set of relative momenta 
qx, so that fll.l2| ) holds. 
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However, there is no one-to-one relation between the momenta p^ and the relative momenta 
q\. We are free to add an extra current, which only flows along the boundary of one of the 
polygons, but never enters or leaves the graph at any vertex. This has no influence on the particles. 
More explicitly, let Xa £ sl(2) be a Minkowski vector for each polygon A E II, and consider 
the following transformation of the relative momenta, 

Q\ i-> <7a + Xa_ a -Xa v (1.13) 



Inserting this into (1.12), we find 



l-> Ptt + E (XA_ A - XA A ) =Pn + J2 XA_ X , ~J2XA x = Pn- (1.14) 

For the first equality, we split the sum into two, and then we replaced the index A G r n in the 
first sum by its successor A' G r w . Thus we shifted the summation index cyclically. But now, 
remember that for two successive links A, A' E r n , we always have A A = A_v- Therefore 
the two sums are equal and cancel. We find that the momentum vectors are invariant under any 
transformation of the form ( |1.13| ). 

Before considering this in more detail, let us show how the link variables can be used to solve 
the constraint equation P = Mj . The total momentum vector can be written as 

p = E p* = E E ?A = E 9a + E <?a + E Qv = - E ( 1 • 1 5 ) 

t 7T ag-Ttt Aer+ Aer_ ^er.oo rjer^ 

The first equality is the definition of P, and the second follows from ( |1.12| ). To derive the third 



equality, we used ( |1.8[ ) and (|1 .9[), which implies 

U r v = r+ u r_ u r^. (i.i6) 



Finally, the last equality in fl!.15| ) follows from the fact that q_\ = —q x for all links \ E T. The 



sums over r + and T_ cancel, and the sum over can be written as a sum over i"^. Now, 
consider the equation P = M70. This is obviously satisfied if we choose the relative momentum 
vectors for the external links to be proportional to 70, hence 

q ri = -M v j , q- r) = M v y , r] E T^. (1.17) 

The variable M v is called the energy of the external link r q E r^. The total energy is then the 
sum of all energies assigned to the external links, 

P = £M„ 70 M X M lr (1.18) 

So, what we found is a partly redundant but complete parameterization of the momentum vectors 
p w of the particles, which are subject to the equation P = M70. The independent parameters 
are the relative momentum vectors q\ for A E and the energies for 77 G i^o. 



The redundancy transformations are still given by ( |1.13| ). But there is now a certain restriction 



on the parameters xa- The external relative momenta q v for 77 G must be proportional to 
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7o> before and after the transformation. This is the case if and only if xa = — ^a7o for all 
non-compact polygons A G IJ^, where cu A £ This implies that the energies transform as 

M v ^ M v + U A _ V -C0 Av . (1.19) 

As a cross check, let us count the number of independent redundancies. There are three com- 
ponents of the vector xa for each compact polygon A G 77 , and one real number to A for each 
non-compact polygon A G TJqo. Thus all together we have 3p + Poo parameters. However, 
there is actually one degree of freedom less. If we set xa = — ^7o for all polygons A G Ft, and 
to a = uj for all non-compact polygons A G Ft^, where w G lis some fixed number, then the 
transformations ( |1.13| ) and fll.l9j ) are void. Therefore, the number of independent redundancies 
is 

3p + Poo - 1 = (3£ + 4o) - (3n - 2). (1.20) 

The equality follows from ( |1.1 1| ). The result is what we have to expect. On the right hand side, we 
have the difference between the number of independent link variables, and the original number 
of momentum variables. There are 3£q independent components of the relative momenta q\ for 
A G r + , and £oo independent energies M^ for rj G .Too. Originally, we had 3n independent com- 
ponents of the momentum vectors p w , and we imposed two constraints, the spatial components 
of the equations P = Mj . 



Relative positions 

The positions x n of the particles are subject to the constraint J = S"yo- They can now be 
parameterized in a similar, but in a sense dual way. For each link A G r, we introduce yet 
another vectors z x G sl(2). For an internal link, it represents the relative position of the two 
particles connected by the link, 

z\ = x nx - x^_ x z_ A = -z x , A G r . (1.21) 

For an external link, it is a spacelike unit vector, which specifies the direction of the link in the 
embedding Minkowski space. Remember that we required the external links to be spatial half 
lines. The vector z v for r] G is therefore a spacelike unit vector which is orthogonal to the 
7o-axis, and the same applies to the vector z^, which points into the opposite direction. Using 



the definition flA.6Q we can write 

z v = 7(0i?) » z -v = -7(&j)> V G ^oo, (1-22) 

where (fi v is the angular direction of the external link r\. In figure [Z] we have cut the ADM surface 
apart, and embedded the individual polygons into some auxiliary Minkowski space, to visualize 
the definition of the vectors z\ for internal and external links. 

Finally, there is yet another link variable that can be assigned to every external link 77 G F^. 
Consider a clock which is located at the far end of the link, and which shows that absolute time 
in the centre of mass frame. The time shown on this clock is then equal to the 70-coordinate of 
the particle 7r_^ to which the link is attached, thus 



T v = x° , r] G -Too. (1.23) 
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Figure 2: The geometry of the individual polygons is determined, up to smooth deformations, 
by the relative position vectors z\ for A G Jo, and the spatial unit vectors z v for rj E i~±oo. 
The ADM surface in figure [I] is reconstructed by gluing the polygons together along the 
edges. The relation z_a = — z\ ensures that they always fit together. 



So, we have now assigned a vector z x to each internal link A G r , and two scalars <p v and T v 
to each external link. But these variables are not all independent. It follows from the definition 
fll.2ip that, for every compact polygon, the relative position vectors representing its edges add 
up to zero, 

$>a = 0, AgT7 . (1.24) 

Aer A 

For non-compact polygons, there is a similar relation. Let i], r/ G -Too be two successive external 
links, and let A. = A. v = A_ v > be the non-compact polygon in between. There is then the 
following relation between the clocks and the relative position vectors, which follows from Q1.21[ ) 
and (JQ3D, 

T^-T V + ^4 = Q, AgZToo. (1.25) 

\er A nr 

Hence, the difference between the absolute time on the two external edges is equal to the minus 
the sum of the the time components of the internal edges. 

The statement is now the following. Given the vectors z\ for all A G r, and the clocks T n 
for all r] G -Too, and provided that the consistency conditions fll.24| ) and ( J1.25D are satisfied, then 
the positions x n of the particles are uniquely determined by the system of linear equations fll.21| ) 
and ( |1.23| ), and the constraint equation J = Sj . The proof is not very complicated. Given 
the vectors Z\ for all A G r , can already derive the relative position of any two particles. We 
just have to walk along the internal links from one particle to the other, and add up the relative 
position vectors. The relation ( |1.24| ) ensures that the result is independent of the path that we 
took. 

So, we already know the absolute positions up to an overall translation where 
w G sl(2) is some unknown Minkowski vector. Now, consider the equation J = Sj . Let us 
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express the angular momentum vector J in terms of the link variables. Using the definition fll.6D, 
and inserting the representation ( |1.12| ) of the momentum vectors, we get 



J = lE^'^l = IE E^^aI- ( L26 ) 



Here we used that n = n\ for all A £ r w . It follows from (1.16) that this can be written as 



J = \ E [9a, JCttJ + | E [«a, a^J + § E fe, ^J- (!- 27 ) 

Aer+ Aer_ r/er^ 

The first two sums can be combined, using q A = — q x , and then we can insert the definition 
( |1.21D of Z\. To simplify the last sum, we use the definition fll.171 ) of M v . The result is 



J = \ ^ [<Za, *a] + I E M . [7o, ^_J- (1-28) 
Aer+ agj 1 ^ 

Now, consider an overall translation x w t-^ x n + v. The first term is obviously invariant, and so 
the angular momentum vector transforms as 

j ^ j + lM[ 7o ,«]. (1.29) 

This implies that the spatial components of the unknown vector t> are determined by the equation 
J = S~f . But then we still have the freedom to perform a time translation x n i— > + f7 . To 
fix the remaining degree of freedom, we need to know the absolute time coordinate of at least 
one particle. This additional information is provided by the clocks T^. It is sufficient to know 
one of them, but it does not matter which one. The consistency condition ( |1.25| ) ensures that we 
always get the same results. 

So, all together we find that the link variables provide sufficient information to reconstruct the 
original position and momentum variables x n and p n of the particles, which are subject to the 
constraints P = Mj and J = Sj . The charges M and S can then also be expressed in terms 
of the link variables. According to fll.!8| ) and ( |1.28[ ) we have 



M = ^M„ S = Y,L X , where L\ = |Tr( [q\, Z\] 7 ). (1.30) 

r?eroo \er+ 

The total energy is, in a sense, distributed over the external links, and the total angular momen- 
tum is distributed over the internal links. Note that L_ x = L x , which implies that the sum is 
independent of the decomposition r = r + U r_. 

But still, the parameterization is not unique. We still have the redundancy transformation 
( |1.13p and fll.l9[ ), acting on the relative momenta and energies. Moreover, we never used the 
unit vectors z± rj for 77 e when we derived the absolute positions x w of the particles. We are 
therefore free to perform the following additional redundancy transformations, 

+ e v , V e Too, (1.31) 

where e v e R are free parameters. This is the freedom that we also have when we introduce the 
triangulation in figure |T]. The directions of the external links are not fixed, so within a certain 
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range we can rotated them. This requires a smooth deformation of the ADM surface, but it does 
not affect the positions or the momenta of the particles. 

Taking this into account, we have the following total number of redundancies, 

3po + Poo + 4o - 1- (1-32) 

On the other hand, we also have the consistency conditions to be satisfied by the link variables. 
There are 3p components of the vector equations fll.24[ ), and p^ scalar equations ( |1.25[ ). But 



these are not all independent. Let us add the 7 -components of all equations fll.24] ), and all 



equations ( |1.25| ). Then the clocks drop out, and what remains is the sum of the 7 -components 



of all vectors z\ for A £ r . But this sum is zero, because for each link A £ i~o, there is also a 
reversed link —A £ r . So, the total number of consistency conditions is 3p + p^ — 1. 

The idea is to consider these consistency conditions as first class constraints, and the redun- 
dancy transformations as the associated gauge symmetries. To treat the rotations ( |1.31| ) of the 
external links in the same way, we shall add an auxiliary angular momentum L n for r\ £ 
to the link variables, which is canonically conjugate to the direction </> r/ , and impose the con- 
straint L v = 0, which then generates the gauge transformation Q1.31| ). With this extension of the 
phase space variables, the number of consistency conditions, or constraints, becomes equal to 
the number of redundancies, or gauge symmetries ( |1.32[ ). 

The extended phase space 

Let us summarize all this in the following definition of an extended phase space Q, for a system 
of n free point particles in the centre of mass frame. Since there are finitely many possible ways 
to define a graph r, the extended phase space Q consists of a finite number of disconnected 
components Qp. There is one component for each possible graph r with n particles. We shall 
discuss these global features of the phase space Q in a moment. Let us first consider a fixed 
graph r, and define the component Q r . It is spanned by the following link variables, 

- a relative momentum vector q\ £ sl(2) and a relative position vector z\ £ sl(2) for every 
internal link A £ r , and 

- an energy M v , a clock T, q , an angular momentum L n , and a direction <f) v for every external 
link. 

The vectors are not all independent. For the internal links we have 

q A = -<?a, z-x = -z\, A £ r. (1.33) 
The independent phase space variables are thus q\ and z\ for A £ r + , where Jo = 7 + U T_ can 



be any possible decomposition. The relation (JL33j) is also valid for external links, if we define 
the associated vectors to be 

q v = -Q-v = ~ M v 7o, ^ = -z_, = 7(^), r] e T^. (1.34) 
The symplectic structure on Q is defined by the following canonical symplectic potential, 

6 = J2(T V dM n + L v d<f> n ) - \ Tr(dq A z x ). (1.35) 

r?eroo \er+ 
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Of course, there is no particular motivation for this at the moment. And in fact, we must show 
that this coincides with the original symplectic structure (JTTT|) . Otherwise the definition does 
not make sense. We still want to consider the same dynamical system. The proof will be given 
further below, so let us for the moment assume that the definition is correct. We can then read 
off the following non-vanishing Poisson brackets. For each external link, we have two canonical 
pairs 

{T„, M„} = 1, {L rp( f> v } = l, ^foo. (1.36) 

For each internal link A G r , the vectors q\ and Z\ are canonically conjugate to each other. 
However, we have to take into account the relations ( |1.33| ), stating that the vectors for A and — A 
are not independent. What we find is 

{qlz b x } = V ab =► {q a . x ,z b x } = { q lz b _ x } = - V ab , \er . (1.37) 

All other brackets are zero, in particular those between variables assigned to different links. We 



should also note that the expression ( |1.35D for the symplectic potential is independent of the 
decomposition T = r + U r_. We have to sum over all internal links with some preferred 
orientation, but it does not matter how this orientation is chosen. If we replace A with —A, both 
q x and z x change their sign, and thus the corresponding contribution to is invariant. The same 



holds for the Poisson brackets (1.37). 



Kinematical constraints 

How is this extended phase space Q related to the original phase space of the free particles, 
restricted to the centre of mass frame? We have seen that the link variables are subject to some 
consistency conditions. The idea is to impose them as kinematical constraints on the extended 



phase space Q. Thus, for each compact polygon, there is a vector valued constraint ( |1.24| ) 



z A = J2*x~o, Aen . (1.38) 



AeA 



For each non-compact polygon, there is a scalar constraint (|1.25|) 



2a = 3V-T„ + X>a»0> AeTIoo. (1.39) 

\er A nr 

And finally, for each external link we impose a constraint which eliminates the artificially intro- 
duced angular momentum L v , 

Jr, = L v ^0, r/G/^. (1.40) 

The subspace Sp C Qr defined by these constraints is called the kinematical phase space. 

It is then straightforward to check that the kinematical constraints are first class, and that the 
associated gauge symmetries are the previously considered redundancy transformations. It is 
convenient first to define a general linear combination. We introduce a set of Lagrange mul- 
tipliers, a vector xa £ s 1(2) for each compact polygon A G II , a scalar G R for each 
non-compact polygon A G 17^, and another scalar e v G R for each external link 77 G r^. Then 
we define 

K = |E Tr (XA^A) + E ZA + E e V Jv O- 41 ) 

Aei7 Aeifoo v^rcc 
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Using the definition of the constraints, and rearranging the summation, this can also be written 
as 

K = \ E Tr( XAA z A ) + £ («a_, - %) T„ + E c, L,. (1.42) 

Here we set %a = — ^a7o f° r A £ so that the first sum also contains the contributions from 
the non-compact polygons. It is then very easy to derive the following non-vanishing Poisson 
brackets, 

{/c,q A } = xa_ a - xa a , a g r , 

{/C, MJ = cu A _, ; - wa,, {£, <M = ?7 e r^. (1.43) 

These are obviously the infinitesimal generators of the redundancy transformations Q1.13| ), ( |1.19| ), 
and ( |1.31| ). We can also see that the constraints are not all independent. Setting e r/ = 0, to a = ^> 



and Xa = 7o f° r some fixed u G R, we find that the linear combination ( |1.42| ) vanishes 
identically, again because the sum over all vectors z x for A G T is zero. We have the following 
relation between the constraints, 

E^a + E^a = 0. (1.44) 

Aeii Aen x 

Everything fits together consistently. We can also recover the original phase space spanned by 
the absolute positions x n and momentum vectors p w of the particles. It is the quotient space 
Sr/~, where two states $2 £ Sr are equivalent, $1 ~ <P 2 , if they can be transformed into 
each other by a redundancy transformation generated by the constraints. 

We can now also check that the symplectic potential (|1.35|) coincide with the original one on 



the kinematical subspace Sr- We start from the original expression fll.ip, and insert the definition 



( jl.l2| ) of the momentum vectors p^, 



= IE MP^x n ) = |E E Tr(q A d^J. (1.45) 

Again, we can rearrange the summation using the identity fl!.16| ), which gives 

O = I E Tr (9A davj + \ E Tr(g A daj w J + \ E Trfa dx^). (1.46) 

Aer+ AG-T_ ^er.^ 

Then we replace the index A in the second sum by —A, so that this also becomes a sum over r + . 
And finally we insert the definition of the relative position vectors Z\, the energies M v , and the 
clocks T v . The result is 

6> = §E Tr(<?A dz x ) - ]T M v dT v . (1.47) 

Ae/+ ^e^oo 

Up to a total derivative, which can always be added to the symplectic potential, this is equal to 
( |1.35P on the kinematical constraint surface Sp C Qr, because there we have L v = 0. 
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Mass shell constraints 



In addition to the kinematical constraints, we also have to impose the dynamical, or mass shell 
constraints, providing the actual time evolution with respect to the ADM time t. They are still 



given by fll.2|), an d the Hamiltonian is a general linear combination thereof, 



n = £ Ck, C n = ±Tr( P7r 2 ) + im, 2 « 0. 



(1.48) 



The subspace Vp C Sp defined by the mass shell constraints is called the physical phase space. 

To derive the time evolution equations, we have to express the momenta p n as functions of 
the relative momenta q x for A E r , and the energies M v for 77 E r^. Consequently, the mass 
shell constraints have non-vanishing brackets with the relative position vectors z x for A E r , 
and the clocks T n for 77 E r^. A straightforward calculation yields the following brackets of a 
momentum vector p n with a relative position vector z x , 



if 



ab 



if a e r n , 
if A e r_ n , 

otherwise. 



(1.49) 



And for the clocks we finds 



aO 



K, T v } = K, T„} = - E {M K , T v } V 

ner-n ree-Tvrnr-oo 

Finally, the brackets with the mass shell constraints are 

p w if 7T = 7T A , 
{C^, 2 A } = \ -p w if 7T = 7T_ A , {C,r, T^} 

otherwise, 



,o0 



if 77 e r_^, 
otherwise. 



(1.50) 



Pi 



if 7T 



7T_ 



otherwise. 



(1.51) 



Using all this, we can easily derive the following brackets with H, providing the time evolution 
of the link variables with respect to the ADM time t, 

z x = {H,z x } = ( nxPnx - ( n _ x p n _ x , % = {H,T V } = C W _,P°_,- ( L52 ) 

This is exactly how we expect the relative position vectors z x and the clocks T v to behave when 
the particles are moving along their world lines. We should also note that this is consistent with 
the definitions ( |1.2ip of z x and Q1.23P of T v , and the original time evolution equations (|L4"|). 



Symmetries 

Let us also consider the rigid symmetries of the system at the level of the extended phase space 
Q. We have seen that the charges are given by 



M 



S 



AG/ + 



(1.53) 



where M v are the energies assigned to the external links, and L x are the angular momenta of the 



internal links, which were given by ( |1.30p . The rigid symmetry generated by M is very simple. 
The only non- vanishing brackets are 



{M,T V } 



T 



T„ - AT, 



(1.54) 
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where AT is some real parameter. The symmetry only acts on the clocks, which are all turned 
backwards by the same amount. In figure HJ, this corresponds to a translation of the ADM surface 
backwards in time, or to a forward time translation of the reference frame with respect to the 
particles. The total energy M is thus the generator of the absolute time evolution in the centre 
of mass frame. This must be distinguished from the ADM time evolution with respect to the 
coordinate time t, which is generated by the weakly vanishing, unphysical Hamiltonian TC. 

To find the symmetry generated by S, let us first evaluate the brackets of the link variables 
with L\. A straightforward calculation shows that the only non- vanishing brackets are 

{L x , z± x } = | [70, z± x ], {L\, q ±x } = I [70, q±\], A G r . (1.55) 

This is the generator of a counter clockwise rotation about the 70-axis, which acts on the vec- 
tors associated with the links A and —A. If we sum over all internal links A G r + , then the 
transformation acts in the same way on all internal links, 

{S, z x } = ~[7o,2a], {S,q\} = 5 [70, 9a], A G T . (1.56) 

The resulting transformation is a rotation ( |A.l 1| ) about the 7 -axis, 

qx „ e A<t> ^ 2 q x e" A ^/ 2 , z x ^ e A ^ o/2 z x A G T , (1.57) 

where A0 is the angle of rotation. But now, we have the following problem. Consider the action 
of this transformation on the polygons in figure 0. All internal links are rotated by an angle 
A(p in counter clockwise direction. But the directions of the external links, specified by the unit 
vectors z± v , are unchanged. As a consequence, there is a certain restriction on the parameter 
A(f>, because for too large angles the non-compact polygons are twisted in such a way that they 
can no longer be realized as spacelike surfaces. 

This can be avoided as follows. We redefine the total angular momentum, including also the 
angular momenta L v of the external links, 



Due to the constraints fll.40p , the charges S and J are in fact weakly equal, thus S ~ J. The 
associated rigid symmetries are therefore equal up to a kinematical gauge symmetry. For J, we 
have the additional non-vanishing brackets 

{J,0„} = 1 {J,z ±rj } = |[7o,2 ±7? ], r/GT^. (1.59) 

The second equation follows from the definition ( |1.34D of the unit vectors z± v . Hence, if we 



replace S by J, then the transformation ( |1.57[ ) also applies to the external links. For the relative 
momenta q± ri this is trivial, because they are proportional to 70 and therefore invariant, and for 
the unit vectors z±^ it follows from ( |1.59| ). The charge J generates a rigid rotation of the ADM 
surface with respect to the reference frame, without deforming it. For such a transformation, 
the parameter A0 can be arbitrarily large, and for A0 = 2n we get a full rotation and thus the 
identity. On the extended phase space Q, it is therefore more natural to call J the total angular 
momentum, but for all physical states its value is of course equal to S. 
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Global aspects 



There are some global features of the extended phase space which we neglected so far. For 
example, until now we have fixed the graph r. But it is not possible to use the same graph for 
every possible configuration of the particles. We therefore have to take into account that the 
extended phase space consists of finitely many disconnected components. The same, of course, 
applies to the kinematical and the physical subspace, 

Q = \jQr, S = \JS r , V = \JV r . (1.60) 

r r r 

On the other hand, the original phase space of the particles was a connected manifold. We 
previously argued that this phase space is equal to the quotient space <Sr/~- But this identity 
only holds locally. Each components of the quotient space only covers a certain region of the 
original phase space, namely that region where the graph r defines a proper triangulation of the 
ADM surface. 

So, what we actually have is an atlas of finitely many charts covering the original phase space 
of the particles. An alternative point of view is to consider the transitions between different 
triangulations as large gauge symmetries. The quotient space «S/~, where the large gauge sym- 
metries are also divided out, is then globally equal to the original phase space. We shall look 
at this global structure of the phase space more closely in section 0, where it is schematically 
indicated in figure |8|. 

It is also useful to keep in mind that following hierarchy of phase spaces, which is very similar 
to the Hamiltonian formulation of general relativity, 

Q D S D V 

i I (1-61) 

S/~ V/~. 

The starting point is the extended phase space Q. The kinematical constraints Ql .38| — |T~40| ) de- 
fine the kinematical subspace S, and the mass shell constraints (|1 .48|) define the physical phase 



space V. At each level, the gauge symmetries associated with the respective constraints can be 
divided out. The quotient space «S/~ is the original phase space, where the redundancies are re- 
moved. And the quotient space 7- > /~ is the reduced phase space, which is the set of all classical 
trajectories, that is all physically inequivalent solutions to the equations of motion. 

Finally, there is another somewhat marginal point, which also has to do with the global struc- 
ture of the phase space. The kinematical constraints alone are not sufficient to ensure that the 
geometry of the ADM surface in figure |T] is well defined. For the polygons in figure ^ to be 
well defined, the sequences of edges defined by the vectors z x must be boundaries of spacelike 
surfaces. A necessary condition is of course that all vectors z x are spacelike. But there are more 
consistency conditions. For a compact triangle, for example, the vector product of two succes- 
sive edges must be negative timelike, for the polygon surface to be spacelike with the correct 
orientation. For higher polygons and also for non-compact polygons there are more complicated 
conditions. 

However, they can all be written as inequalities involving the vectors z\, and this is all we need 
to know. They can be treated at the same level as the positive energy conditions, which we have 



20 



to impose together with the mass shell constraints. They pick out an open subset of the actual 
constraint surface, which is the physically accessible subset. The kinematical phase space Sp is 
this subset of the constraint surface, similar to the physical phase space Vr, which is that part of 
the constraint surface where the positive energy conditions are satisfied. This further restriction 
has no influence on the local structure of the phase space, for example number of gauge degrees 
of freedom. But it imposes some restriction on the parameters of the gauge transformations, as 
otherwise we drop out of the kinematical phase space. 

For example, we cannot rotate the external links arbitrarily far, thus there is a restriction on the 
parameters e v in fll.3ip , as otherwise the external polygons are no longer spacelike. This was also 
the reason why the transformation generated by S as a conserved charge is not well defined for all 
parameters. And there is also a restriction on the multipliers ( w for the mass shell constraints. We 
cannot shift a particle along its world line arbitrarily far, as otherwise the polygons attached to 
this particle are no longer spacelike. Of course, this already follows from the original restriction 
(|T3|), so this restriction has already been present before we introduced the link variables and the 
extended phase space. 



2 Spacetime geometry 

The gravitational field of a massive point particle in three dimensional Einstein gravity is a cone 
with a deficit angle of 8nGm, where < m < Mpi/4 is the mass of the particle and G is 
Newton's constant, which is the inverse of the Planck mass Mpi = l/G. In the neighbourhood 
of the world line, a cylindrical coordinate system (i, r, cp) can be introduced, so that the metric 
becomes 

ds 2 = -dt 2 + dr 2 + (1 - AGmf r 2 dip 2 . (2.1) 

This metric is locally flat for r > 0, and there is a conical singularity on the world line at r = 0. 
A spacetime containing n particles can in principle be covered by an atlas of n such conical 
coordinate systems, and the relative motion of the particles can be read off from the respective 
transition functions. 

And alternative and more appropriate way to introduce coordinates is first to consider the 
spacetime between the particles. Let us remove the world lines from the spacetime. It is then a 
locally flat manifold, which is not simply connected. If a spacetime vector is transported, say, 
along a path that winds around a world line in clockwise direction, then we pick up a Lorentz 
rotation, which is called the holonomy of the particle. It does not depend on the precise path, 
and is therefore a quantity that can be assigned to the particle. If the deficit angle of the conical 
singularity is 8nGm, then the holonomy is a rotation by 8nGm about some timelike axis, which 
is parallel to the world line. 

From the holonomy we can read off both the mass and the direction of the motion of the 
particle. It can be regarded as a generalized momentum [[T3|]. So, we have in this picture a well 



defined notion of a momentum. Massless particles can also be included. For a massless particle, 
the world line is lightlike, and the holonomy is a null rotation, whereas a conical coordinate 
system like ( pTT] ) does not exist. It is also possible to introduce relative position coordinates of 
the particles, very similar to the free particle coordinates in the previous section. So let us make 
all this a little bit more explicit. 
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Figure 3: The ADM surface is a conical surface with n tips. It is divided into a collection of 
simply connected polygons. The polygons are embedded into an auxiliary Minkowski space, 
which provides an atlas of local coordinate charts. The edges ±A of two adjacent polygons 
A±\ are mapped onto each other by some Lorentz rotation. This ensures that the polygons 
can be glued together, forming a spacelike surface, which is embedded into a locally flat 
spacetime with n conical singularities. 



Minkowski coordinates 

Suppose that the spacetime is foliated by a family of spacelike slices, labeled by an ADM time 
coordinate t. Then we would like to know, for example, what the relative position of two particles 
is at a moment of time t. Consider the corresponding ADM surface, and some neighbourhood 
of this surface in spacetime. If the world lines are excluded, then this is a spacelike surface 
embedded into a locally fiat spacetime. In this locally fiat spacetime, we can introduce an atlas 
of Minkowski charts. This can be done as follows. First, we triangulate the ADM surface in the 
same way as before. We introduce a graph r, so that the links \ E T are spacetime geodesies, 
and the ADM surface is divided into a collection 77 of simply connected polygons A £ 77. 

For simplicity, let us assume that the topology of space is that of R 2 . Possible generalizations 
will be discussed later. Then, in the neighbourhood of each polygon, we introduce a Minkowski 
coordinate system. This is possible because the spacetime is locally fiat, and the neighbourhood 
of the polygon in spacetime is simply connected. The polygons become spacelike surfaces, 
which are embedded into some auxiliary Minkowski space, as indicated in figure [3|. Locally, 
this embedding Minkowski space can be identified with the spacetime. Each edge A £ 7a 
of a polygon A represents a geodesic connecting two particles. It becomes a straight line in 
Minkowski space. 

To this we can assign a vector z\ £ sl(2), in the same way as before. It represents the relative 
position of the particles tt\ and 7t_a in spacetime. Vice versa, the vectors z\ determine the 
geometry of the individual polygons in figure ^ There are only the following modifications, as 
compared to polygons in figure The vectors z\ now refer to local coordinates on the spacetime, 
and not to a fixed reference frame. We are allowed to perform coordinate transformations. On 
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the embedded polygons in figure |3], a coordinate transformation acts as a Lorentz rotation and a 
translation, and the vectors z x are also Lorentz rotated. 

Moreover, for each pair of edges A and —A, the vectors z x and z_ x still represent the relative 
position of the same two particles. But now they refer to different local coordinates systems. The 
vector z\ is defined in the chart that contains the polygon A x , and the vector z_ x refers to the 
chart containing the polygon A_ A . There is an overlap region of the two charts, which contains 
the links A and —A. The transition function between the two charts must be an isometry of 
Minkowski space, thus a Poincare transformation. Let us only consider the rotational component 
of this transition function. It is represented by an element g x G SL(2) of the Lorentz group. 

In figure |3], this means that the edge —A of the polygon A x is mapped onto the edge A of the 
polygon A A by a Lorentz rotation given by g x . Of course, g_\ is then the inverse transformation. 
We have the following relation between the relative position vectors and the transition functions, 

g x = gx 1 , z-x = -gx z x g x x , A e r. (2.2) 

This relation also applies to the external edges in figure ^|. The spacelike unit vectors z± v define 
the directions of the external edges in the embedding Minkowski space. And there is also a 
transition function associated with these links, defining the relation between two adjacent non- 
compact charts. 

Now, suppose we are given the vectors Z\ and the transition functions g\ for all A G f. The 
geometry of the polygons in figure |3| is then determined up to smooth deformation of the surfaces 
in the interior. The transition functions tell us how these polygons are to be glued together. This 
determines the geometry of the ADM surface, and the way it is embedded into the spacetime. 
So, we conclude that the relative position vectors z\ and the transition functions z\ specify the 
geometry of space at a moment of time, up to some physically redundant gauge symmetries. The 
deformations are in fact gauge symmetries of general relativity, namely spacetime diffeomor- 
phisms acting on the foliation. 

It is useful to think about the polygons in figure |3] as a deformation of the polygons in figure 
with Newton's constant being the deformation parameter. We can take a limit G — > 0, where the 
gravitational interaction is switched off, if we relate the transition functions g x in a certain way 
to the relative momentum vectors q x of the free particles. Let us define 

gx = e 47rG9A = 1 + 4ttG q x + 0(G) 2 , XeT. (2.3) 

Newton's constant has to show up in this relation, because the vector q x has the physical dimen- 
sion of a momentum, whereas g x is dimensionless. The numerical factor will be explained in a 



moment. If we insert this into (|2.2D, and take the limit G — > 0, then in the leading order in G 
we recover the relations ( |1.33| ). Moreover, for all transition functions we have g x — > 1, so that 
the individual Minkowski charts can be unified, providing one global chart. So, in this limit we 
recover the triangulated ADM surface of the free particles, which is globally embedded into a 
flat Minkowski space. 

It is also useful to look at the way the link variables transform under coordinate transformation 
in the Minkowski charts. In each chart, we can act on the coordinates with some Poincare 
transformation. On the corresponding polygon in figure |3|, the transformation acts as a Lorentz 
rotation and a translation. Since none of the link variables refers to the absolute position of the 
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polygons in the embedding Minkowski space, we can again on the Lorentz rotations and ignore 
the translations. So, let h& G SL(2) be the Lorentz rotation acting on the chart containing the 
polygon A. Then, one can easily verify that the link variables transform as 

9x i-> ^a_ a S'a hZl, z\ >-> h Ax z x h A l. (2.4) 

Note that the transition function g\ sees the transformations acting on both polygons Aa and 
A_a, whereas the vector z\ only refers to the polygon Aa and transforms accordingly. 

To see what this coordinate transformation looks like when the gravitational interaction is 
switched off, we replace the parameter h A by a vector xa, so that 

h A = e 47rGxA = l + 4nGxx + 0(G) 2 , A G 77. (2.5) 



In the limit G — > 0, the Lorentz rotation fl2.4fj then reduces to the redundancy transformation 



( |1.13D for the relative momentum vectors of the free particles, again in the leading order in G. 
The coordinate transformations in the Minkowski charts are the redundancies of the interacting 
particle system, analogous to the redundancies of the relative momenta for the free particles. 

Holonomies 

When the polygons are glued together, then at each vertex of the triangulation a conical sin- 
gularity arises. Consider a spacetime vector which is transported once around a particle n, in 
clockwise direction. As indicated in figure ^, the vector is defined in some polygon A adjacent 
to the particle n. The path winds once around the particle. Whenever it hits a link —AG fL^ 
beginning at n, then it continues from the link A G r n ending at n, in the next polygon. The 
sequence of in which the links are crossed is given by the cyclic ordering of the set r w . When 
the vector is transferred from the link —A to the link A, hence from one chart to another, then we 
have to act on it with the transition function g x . 

The holonomy u n G SL(2) of the particle n is the product of the transition functions g x for all 
A G r n . The ordering of the factors is defined by the cyclic ordering of the set r n . Which factor 
is the first depends on the polygon A in which the path begins and ends. Let us introduce the 
following notation. If ir is a particle, and A is a polygon adjacent to it, then T^a = {A', . . . , A} 
is the ordered set of all links ending n, so that the first and last elements A, A' G r n are those 
successive links ending at n, which enclose the polygon A = A A = A_y. We say that the cyclic 
ordering of the set r w is broken at the polygon A. 

With this notation, the holonomy of the particle n, evaluated in the polygon A, becomes 

Wtt.a = Yl 9x (2.6) 

There are as many holonomies u^a of the particle n as there are polygons A adjacent to the ver- 
tex 7r. They are, however, just representatives of the same physical object in different coordinate 
charts. To see this, consider the behaviour of the holonomy under Lorentz rotations ( |2.4j ) of the 
Minkowski charts, 

Utt.A ^ 1*7r,A ^A 1 - (2.7) 

Hence, u^a transforms according to the Lorentz rotation acting on the polygon A. Moreover, 
consider the representatives of the holonomy u^^ x and u t ,a_ a m the two polygons, sharing a 
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Figure 4: There is a conical singularity at each vertex of the triangulation. Its holonomy u n , 
evaluated in the polygon A, is the product of the transitions functions g\ for A G A-,A- m 
the given example, we have r nA = {Ai, A2, A3, A4, A5}, because — Ai and A5 are the edges 
of the polygon A. 



common boundary A G r v and —AG r_ n . For these two representatives, the definition ( fZq ) 
differs by the position of the factor g x . For A = A A it is the last factor, whereas for A = A_ A it 
is the first factor in the product. This implies 

■"tt.Aa = 9x l u nA _ x 9\ (2.8) 

But this is just the statement that u 1T ^_ x and ■u 7r ,A A represent the same physical object in two 
different coordinate charts, where g\ is the transition function between the two charts. Note that 
if there is only one link A attached to the particle tt, then there is also only one polygon A, and 
the relation Q2.8D is just a trivial identity, because u nA = 9\ and A = A a = A_a. 

According to the arguments given in the beginning, we can now read off the mass of the particle 
and the direction of motion from the holonomy. The direction of motion is the axis of the Lorentz 
rotation represented by the holonomy. This axis can be defined by projecting the group element 



u n onto a vector p n , which we call the momentum vector of the particle. According to ( |A.12| ), 
the projection SL(2) — > sl(2) is defined by expanding a group element in terms of the unit 
and the gamma matrix, and then dropping the term proportional to the unit matrix. To give the 
momentum vector the correct physical dimension, we introduce Newton's constant once again, 
and define 

u nA = u n l + A-kGpI a 7 oj p nA = pl A 7 a . (2.9) 
By definition, the momentum vector commutes, as a matrix, with the holonomy, and therefore 



defines the axis of the Lorentz rotation in Minkowski space. It transforms in the same way ( p77| ) 
under coordinate transformations, and between the representatives in the different charts we also 
have the relation (^8|), with replaced by p n . 



In the limit G — ► 0, and with ( |23| ) inserted, the product in ( ^6| ) can be replaced by a sum. 



In the leading order in G, we recover the definition Q1.12D of the momentum vector of the free 
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particle, which is the sum of the relative momentum vectors. To explain the numerical factors 
appearing together with Newton's constant, let us derive the relation between the holonomy and 
the mass of the particle. For the mass to be m^, the holonomy must be a rotation by SuGm^ in 
clockwise direction about some timelike axis. All such element of the group SL(2) belong to the 
conjugacy class of e 47rGm ' r " ) ' , which represents a clockwise rotation about the 7 -axis. 

This conjugacy class of SL(2) is specified by the following mass shell and positive energy 
condition. The scalar u n in (|2Ty|), which is half the trace of the holonomy, must be equal to 
cos(ATrGm w ), and the 70-component of the momentum p n must be positive to fix the direction 
of the rotation. Hence, 



cos(4nGm T ), p° = ±Tr(u^7°) > 0. (2.10) 



Note that the trace of the holonomy is a scalar. It is the same in every polygon, because the trace 
is invariant under cyclic permutations of the factors, 



±TrK) = |Tr(n 9x), (2.11) 



2 



and it transform trivially under coordinate transformations. The mass shell condition is therefore 
a well defined, coordinate independent equation. Actually, there are some subtleties regarding 
the global structure of the group manifold SL(2). But these can be ignored here, and we can take 
( JZ.IOD as a condition for the mass of the particle to be m^. A comprehensive discussion can be 



found in | JT3| 1 . 

It is now also possible to include massless particles. For m n = 0, the holonomy becomes a 
null rotation and the momentum vector is a lightlike vector. Thus a massless particle moves with 
the velocity of light. On the other hand, there is an upper bound on the mass, because the deficit 
angle of a conical singularity must be smaller than 2n. The upper bound is Mpi/4 = 1/4G. 
This is also the value where the cosine takes its minimum —1, and the holonomy becomes a full 
rotation by 2n. The allowed range of the mass parameters is 

< < M P1 /4 = 1/4G. (2.12) 

And finally, to explain the numerical factors, let us consider the limit G — > once again. Between 
the trace of the holonomy and the momentum vector we have the relation QA.13[), 



u 2 = 8tt 2 G 2 Tr(p 2 ) + 1. (2.13) 
Therefore, the mass shell condition implies 

1™ / 9n sin 2 (47rGm fe ) 

|Tr(P, 2 ) = - I^G* (2 ' 14) 
In the limit G — > 0, this obviously becomes the free particle mass shell condition (|1.3|). 

The conical infinity 

There is one feature of the free particle system which has so far no counterpart for the inter- 
acting system. This is the reference frame, defined by the embedding Minkowski space. The 
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Minkowski space in figure [3| is now only an auxiliary space, providing local coordinates on the 
spacetime manifold. Since we are free to perform coordinate transformations, rotating the indi- 
vidual polygons, this Minkowski space does not provide a well defined reference frame. On the 
other hand, a reference frame is needed if we want to set up a proper Hamiltonian formulation in 
general relativity [|TP[]. 

To define a reference frame, we have to look at the asymptotic structure of the spacetime at 
spatial infinity. Let us assume that, far away from the particles, the spacetime looks like the 
gravitational field of a single particle. We'll see in section |], that this can be formulated as a kind 
of asymptotical flatness condition, hence a fall off condition imposed on the metric at infinity. If 
this is the case, then an observer at infinity effectively sees a fictitious centre of mass particle, and 
the rest frame of this particle defines the centre of mass frame of the universe. We may therefore 
identify the reference frame with the centre of mass frame, in the same way as we did this in the 
previous section. 

The gravitational field of the fictitious centre of mass particle is a cone like (j2~7T|). However, 
we have to allow this particle to have a variable mass M, which represents the total energy of 
the universe. And it also receives a spin S, which represents the total angular momentum of 
the universe. It is then possible to introduce a conical coordinate system (T, R, 0). It covers a 
certain region of the spacetime, far away from the particle, so that the metric becomes that of a 
spinning cone [Q], 

ds 2 = -(dT + 4GSd(f)) 2 + dR 2 + (l-4GM) 2 R 2 d<j) 2 . (2.15) 

The spinning cone has a deficit angle of 8tiGM, and a time offset of 8nGS. There is an upper 
bound on M which is the same as ( [2.1 2[ ). For convenience, let us also assume that M is positive. 



This is actually not needed, but one can show that the positive energy conditions for the particles 
imply that also the total energy M is positive. Hence 

< M < Mpi/4 = l/AG. (2.16) 

The conical coordinate system (T, R, (j>) can then be regarded as a fixed reference frame. It 
replaces the Minkowski frame of the free particle system. The radial coordinate R defines the 
distance from the fictitious world line of the centre of mass, the time coordinate T defines the 
absolute time in the reference frame, and the angular coordinate (ft, which has a period of 2%, 
defines the angular orientation of the reference frame. If we take the limit G — » 0, with M and 
S fixed, the spinning cone becomes a flat Minkowski space, and the fictitious world line of the 
centre of mass is the 70 -axis. 

So, the reference frame is also a deformation of its free particle counterpart. There is only 
the following technical difference. The reference frame does not provide a global coordinate 
system, which covers the whole spacetime. It only provides another local chart, which covers a 
neighbourhood of infinity. This is the region outside a cylinder surrounding all the world lines. 
To define the absolute positions of the particles with respect to this reference frame, we cannot 
just read off their coordinates. The particles are not inside the chart where fl2.15| ) applies. We 



have to define them indirectly. At each moment of time, we have to tell how the ADM surface is 
embedded into the spinning cone at spatial infinity. 

Given the geometry of the ADM surface, we have to fix two additional parameters to define 
its embedding into the spinning cone, because there are two Killing symmetries of the metric 
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fl2.15p . These are the time translations T i— > T — AT and the spatial rotations <p i— > + 
Hence we have the same rigid symmetries of the centre of mass frame as before. The symmetries 
are the possible rotations and translations of the reference frame with respect to the rest of the 
universe [|19|]. To fix this freedom, we have to assign some additional variables to the external 
links, telling us how they behave when the ADM surface is embedded into the reference frame. 

Since we are still free to deform the ADM surface smoothly, without affecting the locations 
of the particles, let us impose the following restriction on the external links. As for the free 
particles, we require them to be spatial half lines. A spatial half line in the spinning cone is a 
spacelike geodesic extending to infinity, which is orthogonal to the fictitious axis at R = 0. An 
alternative definition is to say that a spatial half line is orthogonal to the Killing field of time 
translations. This is obviously a straightforward generalization of the definition in Minkowski 
space, where a spatial half line is orthogonal to the 7 -axis. 

One can then easily show that on every spatial half line the conical coordinates T and <p then 
converge in the limit where R goes to infinity. So, on each external link we have 

T^T V , 4>^4> v , R -> oo, T] G -Too. (2.17) 

The physical interpretation of the link variables T v and 4> v is almost the same as before. Let us 
consider an observer sitting at the far end of the link r\, hence at spatial infinity. The variable 
(prj tells us where this observer is, that is in which direction, and the variable T v represents the 
absolute time at the location of this observer, hence it defines a clock. 

On the other hand, given the geometry of the ADM surface, specified by the link variables z\ 
and g\, and also the variables T v and (p v , we know how to embed the space manifold into the 
spinning cone, and thus we have complete information about the state of the interacting particle 
system at a given moment of time. We know the relative positions of the particles with respect 
to each other, and the absolute positions with respect to the reference frame. 



Consistency conditions 

There are also some consistency conditions to be satisfied by the link variables, in analogy to 
the kinematical constraints for the free particles. First of all, for a polygon A in figure |3| to be 
well defined, the sequence of edges given by the vectors z\ for A G Ta must be the boundary of 
spacelike surface. For each compact polygon we must have 

£z A = 0, AeI7 . (2.18) 

Aer A 

Additionally, there are some inequalities to be satisfied, for example all vectors z x must be 
spacelike, and the boundaries must be correctly oriented. The precise form of these inequalities 
is still not important, for the same reasons as before, which we discussed in the very end of 
section |XJ. 

Then there are also some consistency conditions for the ADM surface to fit into the spinning 
cone at spatial infinity. It must be possible to embed the non-compact polygons into the spinning 
cone, so that the edges fit together correctly. We have indicated this, somewhat schematically, 
in figure |[ Each non-compact polygon A G TI^ covers a certain segment of the spinning 
cone, which is bounded by the external links r] and rf '. They are spatial half lines with conical 
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Figure 5: The non-compact polygons are embedded into the spinning cone. The spinning 
cone on the right defines the reference frame, whereas the Minkowski space on the left only 
provides local coordinates. The vectors z± rj are related to the conical directions 4> v of the 
links in the spinning cone, in such a way that the average orientation of the Minkowski frame 
is the same as that of the conical frame, so that effectively the embedding Minkowski space 
can also be regarded as a reference frame. 



coordinates ( [2 . 1 7| ) . For the polygons to fit in, there must be some relation between the vectors 
z\, defining the geometry of the polygons in the embedding Minkowski space, and the external 
links variables T v and (fi v . 

To find this relation, let us impose one further restriction. We are still free to perform arbitrary 
Lorentz rotations in all Minkowski charts, including the non-compact ones. On the other hand, 
each non-compact chart has some overlap region with the conical chart at infinity. Hence, there 
is a transition function between the Minkowski coordinates and the conical coordinates. This 
transition function is a local isometry between the spinning cone and Minkowski space. There 
is a special class of such isometries. They have the property that the axis of the spinning cone is 
mapped onto the 7 -axis of Minkowski space. The most general such local isometry is given by 

(T,R,<f>) i > (T + TA (0)) 7o + fl 7 (0- aA (0)), (2.19) 

where a and r are two linear functions satisfying 

r A (0) = AGS a' A (<t>) = 4GM. (2.20) 

Note that this is only a local isometry, because it does not respect the periodicity of (p. There is 
a two parameter family of such local isometries, because we are free to add constants to r A and 
a a- Clearly, this corresponds to a time translation and a spatial rotation, either in the spinning 
cone or in Minkowski space. 

Let us choose the Minkowski coordinates in the non-compact charts so that the transition 
function to the conical coordinates is given by such a local isometry. Hence, for each non- 
compact polygon A 6 TI^ there exists a pair of linear functions ta and a a, so that in the 
overlap region of the respective Minkowski chart with the neighbourhood of infinity the transition 
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function ( |2.19| ) applies. What can we then say about the unit vectors z± v and the transition 



functions g± v associated with the external links? 

Consider first the transition functions. Each pair of external links ±77 represents a triple over- 
lap region between two Minkowski charts, containing the polygons A ±7? , and the conical chart, 
which provides a globally defined chart in a neighbourhood of infinity. There is thus a certain 
relation which involves the transition functions g± v between the two Minkowski charts, and the 
transition functions fl2.19| ) to the conical coordinates. The Minkowski coordinates in the two 
adjacent charts differ by a time translation, which is given by the difference ta a (0) — ta_ a (0)> 
and a spatial rotation, which is given by the difference ta a (0) — ta_ a (</>). 

Note that these differences are constant, because in both polygons the functions ta(0) and 



«a(0) are linear according to fl2.20Q . The translation can again be ignored. The rotation tells us 



that the transition function between the Minkowski coordinates is 

9, = e- 4lGM ', g. v = e^ GM ^ i,6r„, (2.21) 
where the energy M, n is given by 

8nGM v = a A _ v (0) - a Av (0) . (2.22) 
This can be evaluated at any point on the links ±77, again because of Q2.2(J| ). We can also take the 



limit ( |2.17[ ), where the definition of M v becomes 

8ttGM v = 9+ - 9~. (2.23) 



The abbreviations 9^ are defined as 



-aiA v ((f> v ), 9 V = -a A _ v (<j>r,)- (2.24) 



Of course, in the limit G — > 0, the relation Q2.21| ) reduces to the free particle relation Q1.17| ), 



stating that the external relative position vector q x is proportional to 70 . To see that the variable 
M v can still be interpreted as an energy, consider the definition ( |2.23[ ). It is the usual relation 
between the energy of a particle and the deficit angle, if it is measured in a frame which is not the 
rest frame of the particle [ [13] ] . And in fact, the difference on the right hand side defines a certain 
angle, which can be regarded as a deficit angle. 

To see this, we have to consider the unit vectors z± v , defining the directions of the external 
edges in figure |5|. First of all, we observe that a local isometry ( p.!9p maps a spatial half line in 



the spinning cone onto a spatial half line in Minkowski space. So, the unit vectors z± v are still 
orthogonal to the 70-axis. To find the directions of these vectors, we have to evaluate the right 
hand side of ( |2.19| ) in the limit ( [2.17D , on the external edges ±77. This gives 



Zr, = l{<Pv + C> z -v = -7(0r, + e~ ), Tie 7^. (2.25) 

Obviously, the difference 9+ — 9~ is the angle between the edges — r] and tj in the embedding 
Minkowski space, on the left hand side of figure |5J. We may call this the deficit angle of the 
external link 77. It goes to zero in the limit G — >• 0. 

The variables 9^ are called the deviations. They tell us how much the directions of the vectors 



(2.25) in Minkowski space, on the left hand side of figure pi deviate from the directions of the 
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external links in the spinning cone, on the right hand side of the figure |5|. There is no coun- 
terpart of these variables for the free particles, because the spinning cone as a reference frame 
is then replaced by the embedding Minkowski space itself, and therefore the deviations are all 
zero. However, we shall now show that here the deviations are also not independent. They are 
determined by the other external link variables. 



We already have the relation ( |2.23[ ). Another relation follows immediately from the definition 



< J2.24| ). If we replace r] by rf in the second equation, and use that A r/ = A_^/ for two successive 



external links rj, rj' e i"^, then it follows from ( ]2.20D that 

Q- -e+ = -AGM (</y - <j) v ). (2.26) 

This relation has a simple geometric interpretation in figure |5|. For the polygon A to fit into the 
spinning cone between the spatial half lines with conical coordinates <p ri and </y, the opening 
angle, hence the angle between the external edges i] and —rf in Minkowski space, must be equal 
to 1 — AGM times the opening angle in conical coordinates, thus the difference 0^ — 0,,. This is 



the factor that shows up in ( [2.1 5\ > in front of the angular coordinate, thus it relates a conical angle 
to a metric angle. 

Now, suppose we are given the energies M v and the directions 4> v of the external links in 
the spinning cone. Then we can solve the system of linear equations Q2.23| ) and fl2.26| ) for the 



deviations. We have 21^ equations and 2i 00 unknown variables. However, the equations are not 
all independent. If we add them all, then the deviations drop out. The differences (fi v / — <\> n add 
up to 2ir, and what we get is 

M = J2 M v . (2.27) 

Hence the total energy is again the sum of the energies of the external links. Or, in a geometric 
language, the deficit angle of the spinning cone on the right hand side of figure |5] is the sum of 
the deficit angles of the external links on the left hand side. 

So, given the energies M n and the conical directions 0,,, the deviations are determined only up 
to an overall constant 9^ 1— > 9^ + 9, where 9 E E is an unknown angle. In figure [5], this means 
that we are still free to perform an overall rotation of all non-compact polygons in the embedding 
Minkowski space on the left hand side. We need one extra condition to fix this freedom, so that 
the deviations are given uniquely as functions of the other link variables. And of course, we want 
that in the limit G —> we get 9± — > 0. 

The extra condition that we impose is the following average condition. We want that, in a 
certain sense, the overall orientation of the Minkowski frame on the left hand side in figure |5] is 
the same as the overall orientation of the spinning cone on the right hand side. This condition 
can be formulated as follows. Consider a certain angular direction in the spinning cone, not 
necessary the direction of one of the external links. A radial line pointing into this direction is 
mapped onto a radial line in Minkowski space, which points into the direction — a a (0), where 



A is the polygon that contains this radial line. This can be read off from ( [2 . 1 9| ) . 

Hence, «a(0) tells us how much a particular angular direction in Minkowski space deviates 
from the corresponding angular direction in the spinning cone. The average condition now re- 
quires that on average this deviation is zero. More precisely, the integral Qa(0) over all angular 
directions vanishes. Hence, 

r"'d0a A (0) =0. (2.28) 



31 



Here we took into account that each non-compact polygon A G 11^ covers a certain range of the 
conical coordinates, namely the one between the external edges 77, 7/ G r^, where A = A r/ = 

This average condition imposes one extra restriction on the Minkowski coordinates in the non- 
compact charts. To see that it fixes the values of the deviations, let us simplify it a little bit. Since 



a a is a linear function, we can evaluate the integral in ([2.28D . The value of the integral of a linear 



function is half the sum of the values of function at the end points, multiplied by the length of 



the integration interval. With the irrelevant factor of one half dropped, and the definitions Q2.24D 
inserted, we get 

AG/Too 

We can then also make use of the relations ( J2.26D and Q2.23D , and rearrange the summation 
slightly, which gives the equivalent condition 

SM,(8+ + e-)=0. (2.30) 

We see that this provides an additional linear equation, which fixed the absolute values of the 
deviations. In the limit G — > 0, it is easy to verify that the unique solution to the system of linear 
equations ( H23l ), flESft , and flOOl ) is 6± = 0. 



So, after this somewhat technical derivation, what is the conclusion? The deviations 6^, and 
thus the unit vectors z± rj are uniquely specified by the external link variables (j> v and M v . The 
geometry of the ADM surface and its embedding into the reference frame is finally specified by 
almost same independent link variables as previously for the free particles. We have the relative 
position vectors z\, and the transition functions g\ for A G r + , and the energies M v , the clocks 
T v , and the directions (fi v for A G r^. We shall later also introduce an auxiliary phase space 
variable L v , representing the angular momentum conjugate to </> r/ , but for the time being we do 
not need this. 

There is then, finally, another consistency condition which involves the clocks. It ensures that 
the ADM surface can be embedded into the spinning cone, so that the external links are in fact 
spatial half lines with the given conical time coordinates. Consider again a non-compact polygon 
A G i7oo, and the external links 77, 77' G -T^ with A = A,, = A_,,/. Both edges are spatial half 
lines in Minkowski space. The constant 70-coordinates of these edges are given by 

T v + r A ((f) v ), 2V+ta(0d')> (2.31) 



respectively. This follows again from ( [2.19D , evaluated in the limit Q2.17D . The difference be 



tween these coordinates is called the time offset of the polygon A. It is given by 

Trji -T ri + r A {(f>n') - r A {(f> v ) = -T, n + AGS{cj), n > - <f> n ). (2.32) 

Now, the same difference is given by the sum of the 7 -components of the relative position 
vectors z\ defining the internal edges A G i~A H r of the polygon A. Hence, we have the 
relation 

TV -T v + 4GS(<P V , - </>„) + £ z x = 0. (2.33) 

Aer A nr 
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This is obviously a generalization of ( |1.25| ), to which it reduces in the limit G — » 0. We can 
use it to compute the parameter S of the spinning cone, thus the total angular momentum of the 
universe. We have to sum over all non-compact polygons in ( |2.33[ ). Then the differences between 



the clocks T v > — T v drop out, and the conical angles 4> v > — (f> v add up to 2ti. The result is 

8nGS = - Yl E 4 (2-34) 

Aen^ \er A nr 

On the right hand side we have to sum over all internal edges of all non-compact polygons. To 
simplify this, we may equally well sum over all compact polygons as well. For each compact 



polygon, the sum over its edges is zero, according to ( |2.18| ). But then the sum just goes over all 
internal links, thus 

5 = -F^E^A = E^, where L x = -J- (z° x + z°_ A ). (2.35) 



This looks very similar to ( |1 .30Q . The total angular momentum is distributed over the internal 



links. To recover the free particle expression for L x = L_ x , we use the relation (|Z2j), which tells 
us that 

Lx = Tr(( &A z x g, 1 - z x ) 7 °). (2.36) 

We can then write the transition function g x as an exponential of the relative momentum vector 
q x , and expand this up to the first order in G as in ([231) . I n me leading order in G, we then 
recover ( |1.30p . 

So, we finally see that not only the geometry of the ADM surface and its embedding into the 
reference frame is determined by the link variables, but also the parameters M and S of the 



spinning cone. And we have the same kind of consistency conditions, a vector equation ( [2. IE 
for every compact polygon, and a scalar equation (|2.33|) for every non-compact polygon. 



Redundancies 

In the beginning we already saw that the redundancy transformations of the free particles system 
are replaced by the coordinate transformations, hence the Lorentz rotations in the Minkowski 
charts associated with the polygons. We found that the link variables transform as 



9x 



hA_ x g\h A l, z x ^h Ax z x h A l, (2.37) 



where h A G SL(2) represents the Lorentz rotation acting on the polygon A. There is now a 
restriction on these parameters, because the Minkowski coordinates in the non-compact charts 
are related in a certain way to the conical coordinates. We are only allowed to perform spatial 
rotations, thus for non-compact polygons A G L 1 ^ we must have 



h A = e - 47r ^7o ; A G iloo, (2.38) 
where cu A G K specifies the angle rotation for the polygon A. This is also the restriction that we 



had for the free particles. The definition fl2.38j ) is chosen so that the energies transform in the 
same way as before, thus 

M v ^ M v + u A _ v - u Av . (2.39) 
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Clearly, the link variables T v and <p v are invariant under such a transformation, because we only 
change the local Minkowski coordinates, but not the embedding of the ADM surface into the 
spinning cone. There is then, however, a transformation of the deviations. For the transformation 



( [2.37D to apply also to the unit vector z± v in ( [2.25D , we must have 



6+ - 8itGuj a 



o: 



9 — 8ttGcua. 



(2.40) 



Now, one can easily check that this is consistent with fl2.23| ) and ( |2.26| ), hence with the definition 
of the deviations as functions of the other link variables. But it is not consistent with ( [2.30D . So, 
there is something wrong with this transformation. 

Before we look at this more closely, let us consider a second class of redundancies. As for 
the free particles, we are free to rotate the external links within a certain range, which includes 
a smooth deformation of the ADM surface. However, due to the screw like geometry of the 
spinning cone, a rotation of the external link 77 not only affects the conical angular coordinate 
(frrj, but also the conical time coordinate T v . It is not difficult to see that a spatial half line in the 
spinning cone is always a line of constant T + AGScf). Thus if we fix the start point of a spatial 
half line and rotate it about this point, then its coordinates, in the limit R — > oo, transform as 



+ e 



T„ i— > TL — AGS e 



'/• 



(2.41) 



where e v G R defines the angle of rotation in conical coordinates. To see how such a transfor- 



mation acts on the embedded polygons in Minkowski space, we use the definition ( j2.24[ ) of the 
deviations, and the property ( |2.20[ ) of the functions a a- Note that these functions are now fixed, 
because we only rotate the external link, but we do not change the transition function between 
the conical coordinates and the Minkowski coordinates. It follows that 



el 



9+ ~AGMe v , 



6: 



6~ - AGM e 



(2.42) 



Together with (|2.41D , this implies that the unit vectors z± v are rotated by (1 — AGM)e v about 
the 70-axis in the embedding Minkowski space. Again, we recover the factor 1 — AGM, which 
relates a conical angle to the corresponding angle in Minkowski space. And we also see that 
there is a certain restriction on the parameters e v , because the non-compact polygons are twisted 
when the external edges are rotated, and for too large parameters they are no longer spacelike. 

One can again verify that the transformations ( |2.42| ) are compatible with ( |2.23| ) and ( |2.26[ ), but 
not with ( [2.30D . To understand this, let consider the infinitesimal generator of the most general 
redundancy transformation. To define the generator of a Lorentz rotation acting on a compact 
polygon A G TIq, we introduce a vector xa G sl(2). It is related to the actual parameter of the 
Lorentz rotation by h A = e 4nGxA . For a non-compact polygon A G /Too, we set xa = — ^a7o 
for some uj/\ G R, which is then also compatible with Q2.38[ ). And finally, the generator of a 
rotation of the external links is defined by e v , which is also the parameter in ( |2.41| ). 

For the transition functions and the relative position vectors we find the following variations, 

5g x = AttG(xa_ x 9x - 9\Xa x ), 5z x = AtvG [xa x , z x ], A G T . (2.43) 
This is the infinitesimal generator of ( |2.37[ ). For the external link variables we find 



5M V = uj a - uja 



5T„ 



-AGSe, 



v g r 



(2.44) 
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which is the generator of ( J2.39D and ( |2.41| ), respectively. And finally, for the deviations this 
implies 

56f = —AGM e n - 8nG u Av , 60~ = -AGM e v - 8nG u A _ v . (2.45) 



In the limit G — > 0, these are precisely the transformations (JTT43J). However, there is one crucial 
difference. For the free particles, we found that a certain combination of redundancy transforma- 
tion was trivial. Setting e v = and xa = — ^7o for all polygons, with uo G K, we found that the 
transformation was void. But here we still get a non-trivial transformation, namely 

5g x = AiiGuj [«a,7o], Sz x = A<kGlo [z x , To], S9± = -%-kGuj. (2.46) 

On the left hand side in figure |5], this is a simultaneous rotation of all polygons in the embedding 
Minkowski space, by SitGu in clockwise direction about the 7 -axis. It is this transformation 
that is not compatible with the average condition fl2.3(jp. 



It is now obvious what this means. Since we wanted the overall orientation of the Minkowski 
frame to be the same as that of the conical frame, we are not allowed to perform such an overall 
rotation. There is a restriction on the parameters uj a and e n . The condition is that the average 
condition ( |2.30P must be preserved. Inserting the variation of the energies M v and the deviations 



7) 



from above, one can easily derive the following restriction to be imposed on the parameters 

cua and e. q , 

E (</y - <t>r,) + Yl e v M v = °- ( 2 - 47 ) 

As for the free particles, we have one redundancy less than there are parameters. It is only a 
technical difference that this is now due to a restriction on the parameters, whereas for the free 
particles a certain combination of the parameters drops out. We shall look at this again at the 
phase space level in section |3], and there we shall see that the analogy to the free particle system 
is indeed very close. 



Inserting and removing links 

There is another class of redundancy transformations, which we more or less ignored so far. Of 
course, we can also change the triangulation of the ADM surface, and this should not affect the 
state of the particles. If we consider the time evolution, for example, then we have to do this 
from time to time, since not every configuration of the particles admits same triangulation. For 
example, if we have a graph r where the particles 7r A and 7r_ A are connected by an internal link 
A, then it is in general not possible to use this triangulation for a state where these particles are 
far apart, with many particles in between. 

What typically happens is that a spacelike geodesic connecting these particles intersects with 
the light cones emerging from the other particles in between. This is actually not a feature intro- 
duced by the gravitational interaction. The problem also arises for the free particle system, where 
we just ignored it. So, let us show how to transform from one triangulation to another. We shall 
later consider this as a large gauge transformation, in contrast to the redundancy transformations 
above, which are smoothly generated gauge symmetries. A convenient way to describe a general 
transition between two different triangulations is to decompose it into elementary steps. 

An elementary step is the insertion or the removal of a single link. Inserting a new internal 
link is very simple. Suppose we have a graph r with an associated collection II of polygons, 
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Figure 6: A compact polygon can be split into two compact polygons by inserting a pair 
of links ±A. The new link variables g±x and z±\ are determined by ( 2.48 ) and ( 2.49 ). 
For the reverse transformation, one first has to perform a redundancy transformation, so that 
g±x = 1. The two Minkowski charts can then be unified, the polygons can be glued together, 
and the links ±A can be removed. 



and let A G 77 by a polygon with more than three edges. It can be either a compact or a non- 
compact polygon. For a compact one, the typical situation in sketched in figure Consider 
one of the diagonals of this polygon, and suppose that the polygon surface can be deformed, so 
that the diagonal is contained in the surface. Then it is possible to cut the polygon into two new 
polygons. 

The insertion is defined as follows. The new graph f = TU {A, —A} is obtained by adding a 
pair of internal links to the original graph r. The new collection 77 of polygons is obtained by 
removing A and adding A a and A_a. The edges 7a a and 7a A of the new polygons are defined 
in the obvious way. The links A is also added to the set r wx , and — A is inserted into i\_ A at the 
appropriate place in the cyclic ordering. All other subsets of links r A = 7a and TV = TV are 
unchanged. 

Regarding the new link variables g±\ and z±\, we have to define them in such a way that 
the geometry of the ADM surface is the same as before. Since both polygons are originally 
embedded into the same Minkowski chart, it is clear that the newly arising transition function is 
trivial, hence 

g x = g x = 1. (2.48) 

The relative position vector z\ is determined as follows. At least one of the newly arising poly- 
gons is compact, even if the original polygon A was non-compact. Without loss of generality, let 
us assume that A A is compact. Hence, 7\ A contains the edge A, and a sequence of at least two 
other internal edges. The new vectors z±\ are then given by 

*a = -*-a = -E 2 - ( 2 - 49 ) 

K£f Ax -{\} 



This immediately implies that the consistency condition ([2.1 8|) for the new polygon Aa is satis 



fied. Moreover, if the original polygon A was compact, the the same holds for the new compact 
polygon A_a- This follows from ( [2.49D and the consistency condition ( [2.1 8p for the original 



polygon. Similarly, if the original polygon A was non-compact, the the new polygon A_a is also 
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non-compact, and the consistency condition Q2.33D for the new polygon follows from the same 



consistency condition for the original polygon. 

So, the conclusion is that, whenever it is possible to deform the ADM surface so that it contains 
the diagonal of a polygon, then it is possible to add this as a new pair of internal links to the graph. 



If we define the new link variables according to fl2.48Q and ( [2.49D , then the new triangulation 
defines the same local geometry of the ADM surface, up to smooth deformations, and thus the 
same physical state. For the free particle system, the insertion of a new internal link works in the 
very same way, we just have to replace ( |2.48| ) by q x = q \ = 0. 

Inserting a new external link into a non-compact polygon is equally straightforward. Suppose 
that we have a graph f and = {. . . , 77, 77', . . .} is the set of external links. Suppose further 
that the non-compact polygon A,, = A- V i has more than one internal edge. Provided that the 
polygon surface can be deformed appropriately, we can then introduce a new external link g 
between 77 and 77', so that = {. . . , 77, g, 7/, . . .}. The original polygon = A_ v i splits into 
two new non-compact polygons A ?? = A_ e and A e = A^ v /. 

The new link variables associated with g are determined as follows. In contrast to the insertion 
of a new internal link, there is now some freedom regarding the direction of the new link. The 
new conical direction <\> Q can be any direction in the range 

(j) v < 0, < </y. (2.50) 

It has to lie between the original links 77 and 77', but we are free to rotate it, because there is no 
fixed end point. But once the conical direction <p s is chosen, all other link variables are fixed. For 
the same reason as before, we have to set 

M e = =*► g e = g Q = l, (2.51) 

The transition function between the two new Minkowski charts is trivial. The clock T e is deter- 



mined by the consistency condition fl2.33| ), which has to be satisfied for both new polygons A e 
and A_ e , thus 

T e = T v - 4GS(0, - 0,) - E 4 = Trf + AGS^ - <f> e ) +J2 4 (2-52) 

Aef A _ e nf Aef A() nfo 



The last equality is the consistency condition ( [2.33D for the original polygon A^ = A^/, so we 
may use either expression to define T e . 

Finally, consider the deviations 9^. Since they are determined by system of linear equations, 
involving all external links, it is not immediately obvious whether inserting a new external link 
affects the deviations of the other external links or not. So, consider the system of linear equations 
( J2.23P , ( |2.26[ ), and ( |2.30| ), after the insertion, and let us write down only those equation that 



involve the new links ±g. First of all, we see that the relevant term in the average condition 
( |2.3(jp drops out, because M g = 0. So, at least this equation is unaffected. 

There are then two equations of the form ( |2.26[ ), namely those that define the opening angles 
of the polygons A± e , 

6-, - 9+ = —AGM (0 V - <t> e ), 6- - 6+ = -4GM (0, - 0,), (2.53) 



37 



and one equation of the form ( |2.23| ), namely that for the deficit angle of the new link, 



6+ -0- = 8nGM g = 0. (2.54) 



First of all, we see that once we chose the conical direction <j) e of the new link, these equations 



three given equations above, we get 



can be used to determine the deviations 9^. Moreover, if we then eliminate the variables 0* from 



6', - 9+ = -4GM - 0„) . (2.55) 



But this is exactly the equation ( |2.26[ ) for the original polygon = A^ v , before the insertion. 



Hence, it follows that the other external links are not affected by the insertion. 

So, we can insert internal and external links, whenever it is possible to deform the polygons 
appropriately so that the new link becomes a spacetime geodesic contained in the ADM surface. 
And what about removing links? Suppose that ±A is a pair of internal or external links, so that 
A A 7^ A_ A . Hence the links do not represent a self overlap region of a single coordinate chart. 
Then it is possible to perform a coordinate transformation in one chart, so that the transition 
function g\ = g \ = 1 becomes trivial. Note that this is not possible if A a = A_a, because then 



the transformation Q2.37D is a conjugation, which cannot be used to transform a group element 
into the unit element. 

So, in order to remove a pair of links, we first have to perform a smoothly generated redun- 
dancy transformation, or fix a gauge, so that the associated transition function becomes trivial. 
Then the link can be removed, and the two Minkowski charts can be unified into a single one. 
We can remove as many links as we like, until there is only a single polygon left, and the graph 
is a tree. Vice versa, the maximal number of links is reached when all polygons are triangles. 
And finally, we can transform between any two triangulations using the given elementary steps, 
and the previously considered smoothly generated redundancy transformations, thus the Lorentz 
rotations of the polygons and rotations of the external links. 



3 Phase space geometry 

After this somewhat technical derivation, let us now come to the actual subject of this article, 
the phase space structure of the multi particle model. In the previous section we have seen how 
the geometry of the ADM surface, and its embedding into the reference frame, is specified by 
a triangulation and a set of link variables. It is possible to set up the same hierarchy of phase 
spaces, which we also encountered in section [I] for the free particles, 

Q D S D V 

i I (3.1) 

The extended phase space Q is spanned by the link variables. The kinematical subspace S C Q 
is defined by a set of kinematical constraints. They are the consistency conditions derived in the 
previous section, and the associated gauge symmetries are the redundancy transformation. The 
quotient space «S/~ is basically the set of all possible geometries of the ADM surface, and all 
possible embeddings into the reference frame. 
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In analogy to general relativity, we can think of the extended phase space Q as spanned by the 
usual ADM variables, the spatial metric and its conjugate momentum, the external curvature. The 
kinematical constraints are then the diffeomorphism constraints, the associated gauge symmetries 
are the coordinates transformation on the space manifold, and consequently the quotient space 
is, roughly speaking, the set of all spatial geometries. Following this analogy, we can then regard 
the mass shell constraints as a the Hamiltonian constraints of general relativity. They define the 
physical subspace V C S, and the quotient space "P/~ is the set of all spacetimes, or the set of 
all physically inequivalent solutions to the equations of motion. 

As for the free particles, the phase space consists of a finite number of disconnected compo- 
nents, one for each possible graph r with n particles, 

Q = \jQr, S = \JS r , V = \JV r - (3.2) 

r r r 

To make the quotient spaces <S/~ and connected manifold, we have to include large gauge 
transformations into the definition of the equivalence relation, and these are of course the transi- 
tions between the triangulations, as defined in the very end of the previous section. In the context 
of general relativity they are analogous to the large diffeomorphism. 

The whole phase space structure can be regarded as a deformed version of the free particle 
phase space, again with Newton's constant as a deformation parameter. At any stage of the 
derivation we can get back to the free particles by taking the limit G — > 0. Using this we 
can see how the phase space structures are actually affected when the gravitational interaction is 
switched on. The most interesting deformation is that of the symplectic structure, which has been 
mentioned in the introduction. The Poisson brackets of the components of the relative position 
vectors are no longer zero, which means that when the system is quantized these become non- 
commuting objects. 



The extended phase space 

The definition of the extended phase space Q is almost the same as before. There are finitely 
many disconnected components Qp. Each component is specified by a graph r with n particles, 
and it is spanned by the following link variables, 

- a transition function q\ E SL(2) and a relative position vector z\ E Sl(2) for every internal 
link A E -T , and 

- an energy M v , a clock T v , an angular momentum L v , and a direction <p v for every external 
link. 

The only modification, as compared to the free particles, is that the relative momentum vector q\ 
is replaced by the transition function g\. There is still a relation between the variables associated 
with an internal link A and the reversed link — A. According to ( |2.2| ) we have 

9 \ = gx 1 , z-x = -gxzxgx 1 , A e r . (3.3) 

Thus actually only half of the internal link variables are independent. There is also a transition 
function g v for each external link, which is given by 

9v = e -*^TD, g_ v = e 4 « GM ^°, V E 7^. (3.4) 
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And finally we have also defined a vector z v for each external link, which is a spatial unit vector 
given by 

Zr, = 7(^ + 0+), Z-„ = -7 (<f) v + 9~), V^Tco- (3.5) 

The quantities 9^, which we called the deviations, are implicitly given by the following system 
of linear equations 

9-, - 9+ = -AGM (</y - (/>„), 9+ - 9~ = 8nGM v , (3.6) 
and the average condition 

Y,M v (9+ + 9-) = 0. (3.7) 



This implies that the relation ( |3.3| ) also holds for external links. Provided that some consistency 
conditions are satisfied by the link variables, which we shall impose below as kinematical con- 
straints, it is then possible to construct the polygons in figure ^ glue them together, and embed 
the resulting ADM surface into the spinning cone, as indicated in figure |5|. This spinning cone 
defines the centre of mass frame, and also the reference frame. The geometry of the spinning 
cone depends on two parameters M and S. They represent the total energy and the total angular 
momentum of the universe, and they are given by 

M = Y / M„ S = Y,L X , where L x = (z° x + z°_ A ). (3.8) 

This is the definition of the extended phase space Qr- It is a 6£o + ^oo dimensional manifold. It 
is a deformed version of the free particle phase space in the sense that the vectors q\ E sl(2) are 
replaced by the group elements g\ E SL(2). Using the relation ( |2T3| ) in the limit G — > 0, we see 
how this deformation disappears when the gravitational interaction is switched off. 



Non-commutative coordinates 

Now need to know what the symplectic structure is, and to say something about the dynamics of 
the particles we also need to know the Hamiltonian. Both will be derived in section |] from the 
Einstein Hilbert action. So, let us here only give a definition without any further motivation. The 
symplectic structure is in fact the most natural one that can be defined on the given phase space. 
It also agrees very nicely with the symplectic structure derived for the single particle model in 
[JT3p, of which this multi particle model is a straightforward generalization. 

The external link variables still provide two canonical pairs (M v , T v ) and (L v , for of each 
external link r] E r^. For each internal link A E r , the pair (g x , z x ) can be regarded as an 
element of the cotangent bundle of the group manifold T*SL(2) = SL(2) x sl(2). The most 
natural symplectic structure is therefore the canonical symplectic structure on this cotangent 
bundle, which is fixed up to a constant. This constant has to be proportional to the inverse of 
Newton's constant, in order to provide the correct physical dimension of the symplectic potential, 
which is that of an action. The numerical factor can then be derived from the condition that in 



the limit G->0we have to recover the free particle expression (1.35). 



40 



All together, this implies that the most natural symplectic potential on the extended phase 
space Q r is 

= £( r i dM„ + A, d0„) - J- £ Trigx 1 dg x z x ). (3.9) 

r- n o7T(jr . ^_ n 



To see that it reduces to ( |1.35[ ) in the limit G — > 0, we have write the transition function gr A as an 
exponential ( JO] ) of the relative momentum vector q x , and expand the result up to the first order 
in G. The given expression is a well defined one-form on Q r , because it is independent of the 
chosen decomposition f = f + U T_. Using the relations (|3.3|), we easily see that 

Tr(gr! dg_ x z_ x ) = Tr^ 1 dg x z x ). (3.10) 

So, we find that the free particle symplectic potential has a natural generalization. But we should 
emphasize that this is just an ad hoc definition. The only way to derive this expression is from 
the Einstein Hilbert action, and this is what we are going to do in section Q 

Let us then derive the Poisson brackets. These are of course the most interesting objects 
when we want to learn something about the quantized model, without actually performing the 
quantization. Nothing particular happens for the external link variables, 

{T V ,M V } = 1, {L v ,<f> v } = l, ^efoo. (3.11) 

More interesting are the internal links. Here we can directly generalize the results from the single 
particle system. For each internal link, we have a pair (g x , Z\), with the symplectic potential 
being the same as (4.3) in [|T3|] . We just have to generalize the Poisson brackets accordingly, 
which are given by (4.11) in p"3[], and insert Newton's constant. The vector z x generates on the 
group element g x a multiplication from the right, and z_ x generates a multiplication from the 
left, 

{g x ,z a x } = 4nGg x r, {g x , z a _ x } = -AnG YQx, Ag4 (3.12) 

To avoid confusion, we shall here and in the following use the vector component notation inside 
the Poisson bracket, and never write brackets with more than one matrix entry. Since the brackets 
between different components of the group elements g x are all zero, it is sufficient to expand the 
vectors z x = z x j a . The vector components themselves provide a representation of the Lorentz 
algebra, 

{z%z b x } = MGe ab c z% \er . (3.13) 

The brackets between the components of z x and z_ A are zero, because the generators of left 
and right multiplication in ( p.!2[ ) commute. The brackets ( |3.1 1| - |3"T3| ) are therefore the only 
non- vanishing brackets involving the basic phase space variables. As a cross check, one can 
verify that the brackets are compatible with (J3.3D. We can either consider g x and z x as a pair of 
independent phase space variables, or g x and z_ A . 

The crucial difference to the free particle system is obviously that the components of the rel- 
ative position vectors z x have non-vanishing Poisson brackets with each other. At the classical 
level, this is just a special feature of the symplectic structure. But at the quantum level it implies 
that the components of the relative position vectors of the particles do not commute. The par- 
ticles are effectively moving in a kind of non-commutative spacetime. We can also say that the 
geometry of the space manifold defined by the polygons in figure |3| becomes a non-commutative 
geometry. 
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What this means to a quantized point particle, and the quantum spacetime that this particle 
effectively sees, has been studied in some detail for the single particle model in [[T3|]. Since we are 
here not going to quantize the model, we can only say that some qualitatively similar effects are 
expected for a multi particle system. The technical details are however more involved, because 
for a proper quantization we also have to solve the various constraints defined below. For a two 
particle system, this can still be done explicitly. And in fact, one finds some interesting features 
of the quantized model. For example, it is impossible to localize the particles at a point in space, 
and it is also impossible to bring them closer together than a certain minimal distance, which is 
of the order of the Planck length [Ell]. 



Kinematical constraints 

The kinematical subspace Sr C Qr is defined in the same way as for the free particles. We have 
to impose the consistency conditions derived in the previous section as kinematical constraints. 



For each compact polygon A G 77 , we found the relation ( [2.1 8[ ), stating that the edges form a 
piecewise straight, closed curve in Minkowski space, 

2a = E z a~0, AG I7 . (3.14) 

Aer A 



For each non-compact polygon A G 11^, there was a consistency condition ( |2.33[ ). It relates 



the 70-components of the vectors representing the internal edges A G i~A PI r to the conical 
coordinates of the two external edges rj G Ta fl and — rj' G i~A H -T-oo, 

Za = T v , -T v + AGS{(j) v , - 0„) + E *° « °> A e n °o- (3-15) 

Aer A nr 

Finally, there is also a constraint which defines the value of the angular momentum L v , which 
is introduced as an auxiliary variable to obtain a well defined phase space. The constraint is a 
deformed version of ( |1 .40p , 

J V = L V + AGSM V « 0, rj G r^. (3.16) 

At this point, there is no particular motivation for this relation between the angular momentum 
L v and the energies M^. This is what comes out in section |], from the phase space reduction 
applied to the Einstein Hilbert action. The only possible motivation for this particular constraint 
is that the associated gauge symmetry is a rotation < J2.41[ ) of the external link 77 in the spinning 



cone, and that in the limit G — > it reduces to the free particle constraint fll.40| ). 

This applies to all the kinematical constraints. There is also a relation between them, which is 
the same as ( |1.44| ). If we add the 7 -components of the vector constraints for all A G I7 , 



and the scalar constraints for all A G 11^, then we get 

E Z% + E Z A = SttGS + zl = 0. (3.17) 

Ag77 Aelfcx) AG-To 

The first equality holds because the conical time differences T v > — T r] add up to zero, and the 
conical angles </y — 4> ri add up to 2n. What remains is a sum of the 7 -components of the vectors 
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z\ for all internal links. But this is just the definition of the angular momentum S in fl3.8p. 
So, there is one independent kinematical constraint less then the number of equations given by 
( P-14| - [JTT6T ). We'll do a precise counting later on. 

All kinematical constraint are first class constraints. This is now less trivial because they 
involve non-commuting phase space variables. But a straightforward calculation shows that they 
form a closed algebra. Every bracket between two constraints is again a linear combination of 
constraints. The only interesting bracket is 

{Z A ,Z b A } = 8nGe ab c Z A} A G /T . (3.18) 

For each compact polygon A, the components of the vector constraint Z A provide a represen- 
tation of the Lorentz algebra. This is quite reasonable, since we expect the associated gauge 
symmetries to be the Lorentz rotations of the polygons in the embedding Minkowski space, thus 
the coordinate transformations in the Minkowski charts. 

To derive the associated gauge symmetries, it is again useful to define a general linear com- 
bination of the kinematical constraints. We introduce a vector valued multiplier xa £ 5 K^) 
for each compact polygon A G 77 , a scalar multiplier co>a G R for each non-compact polygon 
A G /Too, and another scalar multiplier e v G R for each external link 77 G i - ^. Then we define a 
linear combination like ( |1.4ip , 



£ = I E MxaZ a ) +J2co a Z a + J2^ ( 3 - 19 ) 

Inserting the constraints, and rearranging the summation, this can be written as 

^ = § E Tt Ua x zx) - E ( w a, - wa_J T v + Y,£v L v + 

AG-To j/GToo r/GToo 

+ AGS( E^a(0,'-0,)+E^ M O- (3 - 20) 

AGi7oo ^GToo 

Again, we defined Xa = ~ <^a7o f° r non-compact polygons A G n^, because the first sum 
also involves contributions from the non-compact polygons. The first line is the free particle 



expression ( |1.42[ ). The second line is the total contribution that is proportional to S, which 



comes from both ( J3.15| ) and Q3.16| ). It disappears in the limit G — > 0. 

Since S is also a function of the link variables, the constraints are no longer linear, and this 
makes things a little bit more complicated. But we can use the following trick. According to 



the identity fl3.17p , there is some redundancy in the definition of the multipliers. If we make the 



replacement to a >— » ^a + ^ f° r A G 77^, and Xa ^ Xa—^Jo for A G 77 , where uo G R is some 
fixed real number, then the linear combination ( |3.19| ) is unchanged. This can be verified explicitly 



111 ( f3~20| ), where we once again recover the definition of S when we make this replacement. We 



can therefore, without loss of generality, choose the multipliers so that 



E wa(<V-0„)+£ e v M v = 0. (3.21) 

AGi7oo ^GToo 



Doing so, the expression for K, simplifies to 



K = \ E Ti (Xa x z x ) - E ( w a, - w A _J T v + J2 e n L ^ (3-22) 

AG-Tq r^G-Too ^e-Too 
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which is just the free particle expression. Note that this is still the most general linear combi- 
nation of the kinematical constraints, although we are no longer free to choose all multipliers 
independently. 

It is then straightforward to derive the following Poisson brackets of the phase space variables 
with /C. For the internal link variables, we have to apply the brackets ( |3.12| ) and ( p.!3| ), and what 
we get is 

{)C,gx} = 4-7rG(xA_ x g\ -</aXa a ), {£, z x } = 4kG [xa x , «a], A G T . (3.23) 



This is obviously the infinitesimal generator ( |2.43[ ) of a Lorentz rotation of the polygons in the 
embedding Minkowski space. For the external link variables, we use the brackets ( [3.1 1[ ) and find 



{)C,M V } = u A _ v - u Av , 
{IC,L V } = -4GS(lu a _ v -uja v ), 



{)C,T v } = -4GSe v , 
{JC,(f) v } = e v , 7] G Too. 



(3.24) 



These are the infinitesimal generators (|2.44[), with the appropriate transformation of the auxiliary 



variables L v added, so that the constraints ( P-16| ) are preserved. 

We conclude that the vector constraint Z A generates a Lorentz rotations of the compact poly- 
gon A G II , the scalar constraint Z A generates a rotation about the 7 -axis of the non-compact 
polygon A G TI^, and the auxiliary constraint J n generates a rotations of the external link 
r] G .Too in the spinning cone. And we also recover the restriction Q2.47| ) on the parameters of 



these gauge symmetries, which is the same as ( |3.21| ). This ensures that the average condition 
( |3.7| ) is preserved, hence the overall orientation of the Minkowski frame and the conical frame in 
figure |[ 



Mass shell constraints 

The definition of the physical phase space Vr is also analogous to the free particle system, but 
technically again a little bit more involved. We have to replace the mass shell constraints by the 
appropriate deformed versions. First we recall the definition ( [2.6P of the holonomy of the particle 
it, evaluated in an adjacent polygon A, 



(3.25) 



and also the definition (|2.9|) of the momentum vector, which is the projection of the holonomy, 



Ux,A = U n 1 + AllGpl A 7 a , PnA = PlAla- 



(3.26) 



They transform under coordinate transformations in the Minkowski charts according to (|2.7|). 
This is now expressed in the following brackets with /C, 

{JC,U^ A } =47rGr[XA,«7r,A], Ptt,a} = AnG [Xa, Ptt.a] • (3.27) 

The physical phase space Vr C Sp is defined as a subset of the kinematical phase space, by 



imposing the mass shell constraints and positive energy conditions ( |2.10D , 

cos(47rGr7v), pi = |Tr(n^7°) > 0. 



(3.28) 
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As for the free particles, the Hamiltonian becomes a linear combination of the mass shell con- 
straints. To obtain the correct limit G — ► 0, we have to rescale the constraints by a certain power 

H = £f.C„ C, = U '-™ ( tr? m *K . (3.29) 

To see that this provides the correct limit G — » 0, one has to use the relation ( |2.13| ) between the 
scalar and the momentum vector p n , and expand the cosine up to second order in G. 

It is also not difficult to see that the mass shell constraints are still first class constraints. The 
trace of the holonomy u w commutes with /C, and thus with all kinematical constraints. And the 
scalars u w also commute with each other, because they only depend on the transition functions 
gx and the energies M^. To derive the time evolution equations, we have to find the brackets of 
u n with the other phase space variables. The only variables that do not commute with the mass 
shell constraints are the relative position vectors Z\ for A G r , and the clocks T v for rj G r^. 

Consider first the relative position vectors. We have to derive the brackets 

{u n ,z a x } = l{Tr([[g K ),z a x }. (3.30) 
This bracket is zero unless A G r n or —A G r w . Consider the case A G r n first. The transition 



function g x is then one of the factors in the product, and according to fl3.12| ) the action of z\ is 
to insert a factor of AnG~y a into the product, behind the factor gx- So, what we have to do is to 
break up the cyclic product behind the factor g\, and insert a gamma matrix. The result is 

K, z a x } = 2vrGTr( 7 a n 9.) = 2vrG Tr( 7 a u wAx ) = 16vr 2 G 2 p a nAx . (3.31) 

KSfV A, 



If we have —A G r n instead, then we have to use the second bracket in (J3.12|). It is then a factor 



of — 4iiG~Y a , which is inserted into the cyclic product in front of the factor g \. In this case, the 
result is 

K,^} = -27rGTr( 7 a n 9 k) = ~^G Tr( 7 a u nAx ) = -16n 2 G 2 fi Ay (3.32) 
All together, and again in matrix notation, we get 

{Ctt, Z X } = I -p nAx if 7T = 7T_ A , (3.33) 

I otherwise. 

This is formally the same as ( |1.51D , and it has the same interpretation. The mass shell constraint 
C n generates the motion of the particle n along its world lines. The direction of this world line 
is specified by the momentum vector p wA in the chart containing the polygon A. For A = A A , 
this is also the chart in which the vector z x is defined. The relative position vector z x sees the 
motion of the particles n\ and 7t_a- 

There is also a non-trivial action of the mass shell constraints on the clocks T v . They have 
non-vanishing brackets with the energies M v , and thus with the holonomies g± v of the external 
links. It follows from fl3.4[ ) that for r] G and — r qEr_ 00 we have 

{g v ,T v } = -4nGg vl °, {g- v ,T v } = 4nG~f°g- v . (3.34) 
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Formally, T v acts on the transition functions g± v as if it was the 7 -component of a fictitious 
relative position vector assigned to —77. Performing the same calculation again, it follows that 

{C w , T v } = { <a„ = Pl A -n ' lf l = *-v> (3.35) 
1 otherwise. 



In contrast to ( p.33[ ), the case -n = tx v does not occur, because the external link 77 has no end point. 



And furthermore, we can here evaluate the momentum vector p n in either of the two adjacent 
polygons or A_^. The transition function between the two associated charts is g v , and this 
is a rotation about the 70-axis. Consequently, the 70-component of the vector p T is the same in 
both adjacent polygons A,, and A.,,. 

Finally, if we ignore the dependence of the various vectors on the polygons, hence on the 
coordinate charts with respect to which they are defined, then the time evolution equations are 
the same as fll.52p for the free particles, 



Z x = {H,Z X } = CttaPtt* - C-xPtt-x, % = {K T V } = Cr-„P°_„- (3-36) 

And clearly, the physical interpretation of these equations is also the same as before. The particle 
7r moves along its world line, by an amount that is specified by the multiplier ( n . Thereby, the 
geometry of the ADM surface changes, and also its embedding into the spinning cone. If all 
multipliers are positive, this finally provides a proper foliation of the spacetime manifold. 



Symmetries 

In the Hamiltonian formulation of general relativity, every isometry of the asymptotic metric at 
spatial infinity should be realized as a rigid symmetry at the phase space level JT^]. Let us check 
this for the time translations T 1— > T — AT and the spatial rotations 1— > <p + of the spinning 
cone fl2.15| ). And let us also derive the associated conserved charges. By definition, these are the 
total energy and the total angular momentum. First we have to find out how the symmetries act 
on the phase space variables. Then we can check whether these are symmetries of the symplectic 
structure, derive the charges, and finally we have to show that the charges commute with all 
constraints. 

A time translation T 1— > T — AT of the spinning cone does not affect the geometry of space. 
The whole ADM surface is just shifted backwards in time, with respect to the reference frame. 
Only the clocks T v refer to the absolute time coordinate T of the spinning cone, so only they 
transform, 

Trj^Tr,- AT. (3.37) 

Since AT is a constant, this is obviously a symmetry of the symplectic structure defined in fl3.9|). 
The symplectic potential changes by a total derivative O 1— ► © — AT dM, but the two-form 
J? = dO is invariant. The associated charge is also easy to find. It is the previously defined total 
energy M, which has vanishing brackets with all phase space variables, except for the clocks, 

M = ^M V {M,TJ = -1. (3.38) 

It is also immediately obvious that M commutes with all kinematical and dynamical constraints, 
thus the total energy is in fact a conserved charge. 
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The same applies to a spatial rotation <fi i— > <p + A<p of the spinning cone. The local geometry 
of space is unchanged, while the whole ADM surface is rotated with respect to the spinning cone 
by an angle A0 in conical coordinates. As only the angular directions 4> v of the external links 
refer to the conical coordinate <p, let us make the following ansatz, 

<f>n ^<p v + A0. (3.39) 

Now, suppose we apply this transformation. Then the deviations 9^ are invariant. They are im- 
plicitly defined in ( |3~^| ) and (|3.7|), and these equations are invariant under the given transforma- 



tion. But if we then look at the definition ( |3.5D of the unit vectors z±\, we find that these vectors 
are also rotated by an angle A0, but now this rotation takes place in the auxiliary Minkowski 
space, where the polygons are embedded. 

The reason can be seen in figure |5|. If we rotate all external links simultaneously with respect 
to the spinning cone on the right, then we also have to rotate the external edges in the embedding 
Minkowski space on the left. This is because we required the overall, or average orientation of 
the two frames to coincide. But then, if we want the local geometry of the ADM surface to be 
unchanged, we also have to rotate the internal edges of all non-compact polygons. Otherwise the 
polygons are twisted, and the ADM surface is deformed. So, to define a proper spatial rotation 
of the particles with respect to the reference frame, we not only have the rotate the external 
links with respect to the spinning cone. We also have to rotate all non-compact polygons in the 
embedding Minkowski space by the same angle A0. 

But then it is actually simpler to rotate all polygons, including the compact ones. This just 
involves another kinematical gauge transformation. In other words, we have to apply the same 
rotation about the 7 -axis to all Minkowski charts, thus rotate all polygons in figure |]by the same 
angle. And finally, we also have to adjust the transition functions. If we rotate the Minkowski 
coordinates in all charts, then we also have to apply this rotation to the transition functions. All 



together, we get the following transformation of the link variables, in addition to ( |3.39| ), 

gx _> e A ^ 2 g x e" A ^/ 2 , Zx ^ e A ^ 2 z x e~ A ^ 2 , A G f. (3.40) 

For internal links, this defines the transformation of the phase space variables z x and g\. For 
external links, it is consistent with the definitions ( p.4| ) and Q3.5|), and the transformation ( |3.39j ). 



Given these transformations, it is straightforward to check that the symplectic structure ( |3.9| ) is 
invariant. Thus we also have a symmetry of the phase space. To find the total angular momentum, 
let us first consider the internal links, and the definition ( Jj.8p of the angular momenta L\. It is 
not difficult to derive the brackets 

{L x ,z±x} = lho,z ±x ], {L x ,g±x} = lho,9±\], A e r . (3.41) 

The phase space function L x generates a counter clockwise rotation of the internal edges ±A in 
the embedding Minkowski space. This was also the defining property of the angular momenta L\ 
in ( |1.55| ), for the free particles. But now we also have to rotate the external links. The generator 



of the rotation ( |3.39[ ) is the angular momentum L v , 

{L V ,4> V } = 1 fjef^. (3.42) 
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To act on all links at the same time, we just have to sum over all internal and external links. 
Hence, the total angular momentum becomes 



j = E l a + E l t ( 3 - 43 ) 

This is formally the same as ( |1.58| ) for the free particles. To see that it is a conserved charge, we 
have to compute the brackets of J with the constraints. It turns out that the only non- vanishing 
brackets are those with the vector constraints for the compact polygons, 

{J, Z A } = | [70, Z A ] « 0, AG n . (3.44) 

This is again proportional to the same constraint, so that J is at least weakly conserved. This is of 
course sufficient to define a conserved charge. We conclude that all symmetries of the spacetime 
are in fact realized as symmetries of the phase space. 

But how is J related to the parameter S of the spinning cone? It is no longer so that J is 
weakly equal to S. Instead, 

j = ^2 l x + J2 L v~ s - AGS J2 M v = ( 1 ~ AGM ) s - ( 3 - 45 ) 

So, the total angular momentum is not the geometric parameter S, which defines the time offset 
of the spinning cone. It is rescaled by a factor 1 — AGM. 

The definition of the total angular momentum is actually somewhat ambiguous. Of course, S 
is also a conserved charge. It is a function of J and M, at least up to terms proportional to the 
constraints. So, there is also a symmetry associated with S. It must be some combination of a 
time translation and a spatial rotation, thus a kind of screw rotation. And in fact, the definition 
of the rotational Killing vector of the spinning cone is also somewhat ambiguous. The Killing 
vector associated with time translations is uniquely defined as that Killing vector which is parallel 
to the fictitious world line of the centre of mass of the universe. Hence the definition of the total 
energy M is unique. 

However, there are two alternative definitions of the rotational Killing vector, and both are 
quite natural. They coincide only if the spinning cone is actually static, thus if the angular 
momentum is zero. One possible definition is to choose the unique Killing vector which has 
closed orbits of affine length 2n. The associated angular momentum is then J, which has the 
property that a symmetry transformation with parameter A0 = 2ir is the identity. But we may 
also define the rotational Killing vector to be the unique Killing vector which is orthogonal to 
the Killing vector of time translations. On a spinning cone, this Killing vector generates a screw 
rotation, and it is this symmetry which is associated with the conserved charge S. 

It is not clear which of the two is the correct angular momentum, S or J. Which one is the 
more appropriate one depends on the context. For example, if we quantize a particle model like 
this, then it is J whose eigenvalues are quantized in steps of h, because the essential property is 
then the existence of closed orbits of affine length 2ir. On the other hand, if we want to describe 
the conical geometry of the spacetime at infinity, then S is more natural, because it is this which 



appears in the metric (2.15) of the spinning cone, and which defines the time offset 8nGS. 
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An interesting consequence of this interplay between quantization and geometry and the 
charges J and S is discussed in the very end of [jZH]. It turns out that, independent of the mat- 
ter content of the universe, only those conical geometries can be realized in a quantized three 
dimensional universe, which fulfill the following Dirac like quantization condition, 

(l-AGM)S „ 

J -eZ. (3.46) 



n 

Since here were are not going to quantize anything, we shall not go into any details and derive 
this quantization condition. But it is more or less obvious that this follows from the fact that J is 
quantized in integer multiples of h. This result is independent of the non-commutative structure 



of spacetime implied by the brackets (3.13), but it also provides a kind of quantized geometry. 



Parity transformations 



Let us also briefly consider a discrete symmetry. Apparently, the Poisson algebra ( [3.13D is not 
invariant under parity transformations, thus under reflections of space, because the Levi Civita 
tensor shows up. But on the other hand we know that the Einstein equations are invariant, and 
thus for every spacetime there exists a reflected spacetime. It is therefore not immediately obvi- 
ous that this symmetry is also realized at the phase space level. And in fact, for the previously 
considered single particle model this was not the case, due to a somewhat ambiguous definition 
of the reference frame [ |I3| ] . 



For a single particle, we had the same Poisson algebra (B.13[), but there was only one vector z 



representing the absolute position of the particle with respect to the reference frame. Now each 
vectors z\ is the relative position vector of two particles, and there are several such vectors. But 
still, the Levi Civita tensor shows up, and we should therefore check whether the problem is still 
present or not. First of all, we have to find out how the various phase space variables actually 
transform under a reflection of space. We have to start at the very beginning, the definition of 
a triangulation. When we perform a reflection of space, we have to replace the graph r by its 
mirror image. 

This is a non-trivial operation, because the cyclic orderings of the subsets r n , i"^ and i~A 
already refer to the orientation of space. For each particle it, we have to reverse the cyclic 
ordering of the set r n of links ending at w, and the same applies to the set of links pointing 
towards infinity. For each polygon A, we not only have to reverse the cyclic ordering of its 
edges in Ja, but also reverse the orientation of the edges. For example, if a polygon A is 
bounded by the edges Ta = {Ai, A2, • • • , A n } before the parity transformation, then after the 
parity transformation we have i~A = {— A„, . . . , — A2, — Ai}. 

Now, consider the reflection of a polygon in the embedding Minkowski space, which is indi- 
cated in figure 0. The edge —A of the polygon A_ A becomes, after the reflection, the edge A of 
the polygon A A . The edge was originally represented by the vector The polygon is then re- 
flected, and additionally the orientation of the boundary is reversed. Hence, the new vector z x is 
obtained from the original vector z_a by a reflection and an overall change of sign. Let us choose 
the reflection in Minkowski space so that 72 1— > —72, hence a reflection at the 70-71 -plane, which 



is given by (A.14) 



Combining this with an overall change of sign, and deriving the appropriate transformation of 
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Figure 7: A parity transformation acts on a polygon in two ways. The polygon is first 
reflected, and then the orientation of the boundary is reversed. If the link — A is an edge of 
the original polygon A_a, which is represented by the vector z_\, then the reversed link 
A becomes an edge of the reflected polygon A\, which is represented by the transformed 
vector z\. 



the transition functions from (|33|) gives 



Z\ l-> 72 «- A 72, 



9\ >-> 72 0- A 72- 



(3.47) 



The external link variables also transform in a non-trivial way. It follows from the definition < ]378| ) 
of the angular momentum that S — » —5. Thus the reflected spinning cone is spinning into the 
opposite direction. The conical coordinates are reflected by the transformation <\> i— > — </>. From 
this we infer that 



'/• 



M n i — ► M 



— L 



'/• 



(3.48) 



The last transformation follows from S i— > — S 1 , and the constraint equation « defined in 
( J3.16| ), which must be satisfied before and after the reflection. The vector constraint ps in 



Q3.14D is obviously also preserved. And for the constraint m defined in ( P-15| ), we have to 
take into account that the external links t] and rj are interchanged, because the ordering of is 
reversed, so that this is also preserved. 

The symplectic potential ( J3.9D can then be shown to be invariant. The external terms T v dM v 



and L v are obviously invariant under the transformation (p.48[). For the internal links, this is 



not immediately obvious. However, we have 7272 = 1, and therefore 



(3.49) 



But this just means that we have to replace r + by T_ in (|3TP|). And since we know that the 

symplectic potential is independent of the decomposition r = J+ U Jl, it follows that it is also 

invariant under parity transformations. 

Finally, we also have to check the mass shell constraints, hence the traces of the holonomies 
of the particles. In the definition ( |2.1 1| ), each factor g\ in the product is replaced by 72<7-a72, 

and the ordering of the factors is reversed. Again, the factors 72 drop out because 7272 = 1. 

So, the group element under the trace is effectively replaced by its inverse. But the trace of the 

inverse of an element of SL(2) is equal to the trace of the element itself. 
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So, we find that all the phase space structures are invariant under parity transformations. But 
what is then the difference to the previously considered single particle system, where this was 
not the case? The difference is that the definition of the relative position vector z x already refers 
to a fixed orientation of the spacetime. It defines the edge of a polygon, where the boundary is 
traversed in counter clockwise direction. This is why the Levi Civita tensor can show up in the 
Poisson algebra, without breaking the parity invariance. 

Open and closed universes 

So far we assumed that the topology of space is IR 2 . But of course, general relativity allows a 
more general topology. The space manifold can be any two dimensional, orientable surface, and 
it can be either compact or non-compact. Let us for example consider a closed universe, where 
the space manifold is a Riemann surface of genus g. In this case, the graph r = T consists of 
internal links only, and all polygons are compact, thus II = IJ . When the Riemann surface is 
triangulated, then we have the following relation between the number of vertices, links, polygons, 
and the genus, 

£ -p = n + 2g-2. (3.50) 

We can use this to calculate the number of physical degree of freedom of a closed universe. 
The extended phase space Qr is spanned by the independent link variables g x G SL(2) and 
Z\ G -6 1(2) for A G r + , hence 

dim(Q r ) = 6£ . (3.51) 

Then we have to impose the kinematical constraints. For a closed universe, we only have the 
vector constraints Za for A G II , and the associated gauge symmetries are the Lorentz rotation 
of the polygons. There is no relation like ( P-17| ), and there is no restriction on the parameters 
of the Lorentz rotations. Thus we have 3po independent kinematical constraints, and also 3po 
independent gauge symmetries, 

dim(S r ) = 6£ - 3p , dim(«S r /~) = 64 - 6p = 6n + 12g - 12. (3.52) 

This is the number of independent kinematical degrees of freedom, thus the phase space di- 
mension before the mass shell constraints are imposed. To obtain the physical phase space, we 
have to impose n additional mass shell constraints, and subtract n additional dynamical gauge 
symmetries. Hence, the number of physical degrees of freedom is 

dim{V r ) = Q£ - 3po - n, dim(P r /~) = 6£ - 6p - 2n = 4n + 12g - 12. (3.53) 

For each particle we have 4 physical degrees of freedom. The particle is moving in a two dimen- 
sional space, thus we have two position and two momentum coordinates. 

Additionally, the closed universe also has some topological degrees of freedom. Each handle 
of the Riemann surface contributes with 12 degrees of freedom. The topological degrees of 
freedom are the holonomies of the fundamental, non-contractible loops in space. There are two 
such loops for each handle, and regarding the holonomy we have to take into account both the 
rotational and the translational components [Q, |J. Hence, for each non-contractible loop we have 
one element of the six dimensional Poincare group, and all together this makes 12g topological 
degrees of freedom, since the genus g counts the number of handles. 
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Finally, there is a consistency condition, which must be satisfied for the universe to be closed. 
This explains the minus twelve in the end. As an example, consider a sphere with g = 0. Then 
we have 4(n — 3) independent physical degrees of freedom. We need at least four particles 
to get a non-trivial dynamical system. In this case, the space at a moment of time looks like 
a tetrahedron, and we have four independent degrees of freedom. The relative motion of two 
particles already determines the motion of the other two particles. For higher genus surfaces, we 
need less particles or even no particles at all to get a non-trivial dynamical system. 

For an open universe, we may also consider more general topologies. The space manifold can 
be represented as a punctured Riemann surface, thus a surface of genus g with k > 1 points 
taken away. Each puncture represents a region that extends to spatial infinity. For each such 
infinity, there is a conical frame, and a set of external links and non-compact polygons covering 
a neighbourhood of this particular infinity. Moreover, each conical frame has its own total energy 
and total angular momentum, and from each infinity an external observer can see the particles 
and also the other observers. The case that we discussed so far, and which reduces to the free 
particle system in the limit G — > 0, is the special open universe with g = and k — 1. The 
closed universes are also special cases, with k = 0. 

Let us count the number of physical degrees of freedom for an open universe. First we have to 
derive the relation between the number of vertices, links, and polygons. It is the same as ( |3.5(J| ), 
but now we have a vertex not only for every particle but also for every infinity. Moreover, we have 
to distinguish between internal and external links. The external links are now those that extend 
to any of the k infinities. But there are still as many external links are there are non-compact 
polygons. Hence, 

i - po = n + k + 2g - 2, = p^. (3.54) 

The phase space dimension is counted in the same way as before. The extended phase space 
Q r is spanned by 6£ independent components of the group elements g x and the vectors Z\ for 
A G -To, and Ai^ real variables M v , L v , T n , (J) v for rj G r^, thus 

dim(Q r ) =64 + 44o. (3.55) 

Then we have to count the number of kinematical constraints. There are 3po components of the 
vector constraints for A G II , p^ scalar constraints for A G n^, and finally scalar 
constraints J„ for n G r^. But now they are not all independent. For each infinity, we have a 



relation like ( |3.17[ ). This implies that the total number of independent kinematical constraints is 
3jp + poo + ioo — k, and this is also the number of independent redundancy transformations. It 
follows that 



dim(Sr) = 6£ - Sp + + k, 
dim(S r /~) = 6£ - Gpo + 2k = 6n + Ylg + 8k - 12. (3.56) 

This is the number of independent kinematical degrees of freedom. To get the physical degrees 
of freedom, we have to impose n additional mass shell constraints, and subtract another n gauge 
symmetries. This gives 

dim(P r ) = 6i - 3p +2£ 00 -n + k, 
dim(7V~) = 6£ - 6p - 2n + 2k = An + 12 g + 8k - 12. (3.57) 
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Once again, we have four degrees of freedom for every particle, and twelve degrees of freedom 
for every handle. And obviously there are eight degrees of freedom associated with every infin- 
ity. This can be understood as follows. There are two charges M and S, defining the conical 
geometry, which can be different at each infinity. The other six degrees of freedom represent the 
relative orientation of the conical frames with respect to each other. 

Consider for example two infinities and the associated conical frames. Think of them as the 
rest frames of two different observers. When the observers talk to each other, they find out 
that they do not agree about the fictitious position of the centre of mass of the universe, and in 
general they are also in relative motion with respect to each other. Hence, every infinity has its 
own centre of mass frame. The relation between two such frames is, at least locally in phase 
space, parameterized by a Poincare transformation. If we fix one reference frame, then for every 
other there are six degrees of freedom, representing the relative position and orientation with 
respect to the fixed one. 

The minus twelve can now be explained by looking at the special case k — 1 and g — 0, which 
reduces to the free particle system in the limit G — > 0. In this case, the dimension of the reduced 
phase space is 4 [n — 1). This is the number of physical degrees of freedom of the relative motion 
of n particles in a two dimensional space. To get a non-trivial dynamical system, we need at 
least two particles. In general, each additional infinity, and also each handle already provides 
sufficiently many topological degrees of freedom, so that we do not need any particles at all to 
obtain a non-trivial dynamical system. These are the usual, say, cosmological toy models, where 
only the topological degrees of freedom are present, and they were the first quantized toy models 
in three dimensional gravity 0] . 

Large gauge symmetries 

Finally, let us also have a look at the global structure of the phase space. We have seen that the 
extended phase space Q consists of finitely many disconnected components Q r , one for each 
possible graph r with n particles. According to ( |3.55| ), they even have different dimensions. 
It is therefore not possible to consider the components Qr of the extended phase space Q as 
overlapping coordinate charts of a single manifold. However, we can still glue them together in 
a certain way, so that the quotient spaces «S/~ and V / '~ become connected manifold. So far, we 
only considered the gauge symmetries generated by the constraints. 

Two states $>i, $2 £ <Sr, belonging to the same component of the kinematical subspace, are 
defined to be physically equivalent, $1 ~ $ 2 , if they can be smoothly deformed into each other 
by a transformation which is generated by the kinematical constraints. This is equivalent to the 
condition that they lie on the same gauge orbit, that is the same null orbit of the symplectic 
structure on Sp- However, there are also states that belong to different components of S, but 
which are nevertheless equivalent because they define the same geometry of space. They are 
then related by a large gauge transformation. In the end of section we have seen how these 
transformations look like. They can be decomposed into elementary steps, namely inserting and 
removing links. 

From the phase space point of view, the typical situation is shown in figure f| There we see 
three components of the extended phase space Q, corresponding to the graphs r^°\ and 
r^ 2 \ Let us assume that both rW and are obtained from by inserting a new pair of 
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Figure 8: The extended phase space Q consists of finitely many disconnected components 
Qr of different dimension. A large gauge symmetry embeds a lower dimensional component 
Sp of the kinematical phase space into the respective higher dimensional component. The 
embedding is given by the insertion of a new link A, as shown in figure H The image of the 
lower dimensional component in the higher dimensional one is the subspace defined by the 
gauge condition g\ = 1. 



links, ±Ai and ±A 2 , respectively. For example, the new links may be two different diagonals of 
the polygon in figure^. The component Q r (o) is the lower dimensional one, and the components 
Q r (i) and Q r <?) have the same higher dimension. The same applies to the kinematical subspaces 
«S r (o), S F (i), and S r (2), which are submanifolds of the respective components of Q. 

The transition from the graph r^°' to r^ 1 ' is a large gauge symmetry, which maps a state 
E S r (o) onto an equivalent state E <S r (i). And a transition from r^°> to r^ 2 ' maps 
the state E S r (o) onto an equivalent state E S r (z). We can think of an embedding 
of the lower dimensional component <S r (o) into the higher dimensional components <S r (i) and 
S r (2). The image of this embedding is a submanifold of the higher dimensional component 
of the kinematical phase space. It is the submanifold which is defined by the gauge condition 
g\j = 1 in S r (i), and g\ 2 = 1 in <S r (2), respectively. 

A general kinematical gauge transformation can be described as follows. As long as we stick 
to a given triangulation, say the gauge symmetries are generated by the kinematical con- 
straints. We follow a gauge orbit, hence a null direction of the symplectic structure on the kine- 
matical subspace S r w. Now, suppose that we want to remove the pair of links ±Ai. Then, 
we first have to fix a gauge where g\ 1 = 1. Thus we have to follow the gauge orbit to some 
state where this gauge condition is satisfied. This state <P^ E S r (i) is the image of some state 
<P^ E S r {o) in the lower dimensional component, and the large gauge transformation takes us 
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there. 

If we want, we can then perform another large gauge transformation, inserting a new link A2. 
Then we have a state in the component <S r ( 2 ) of the kinematical phase space. There we can 
again perform a smoothly generated gauge transformation, following a gauge orbit provided by 
the kinematical constraints, and so on. All components of the phase space are finally connected 
in this way, and consequently the quotient space «S/~ becomes a connected manifold. That it is a 
proper manifold follows from the fact that all components iSr/~ of this quotient space have the 
same dimension. According to ( |3.56[ ), this dimension only depends on the number of particles 
and the topology, but not on the triangulation. 

However, what we still have to show is that <S/~ is in fact a proper symplectic manifold, and 
that the Hamiltonian, or the mass shell constraints are well defined functions thereon. In other 
words, we have to show that the symplectic structure is invariant under large gauge symmetries. 
This is the case if the symplectic potential on the lower dimensional component S r (o) is equal to 
the one that is induced by the embedding into the higher dimensional components <S r (i) and S r (2) . 
And the same has to apply to the holonomies of the particles. It is not difficult to see that this is 
indeed the case. If the inserted link is an internal link A, then the image of the lower dimensional 
component in the higher dimensional one is the subspace defined by the gauge condition g\ = 1. 
If we look at the symplectic potential fl3.9|), we find that on this subspace the terms involving the 
link A are just absent. 

Hence, the symplectic potential is the same as if the link was not there. The same holds if the 
inserted link is an external link 77. In this case, the image of the lower dimensional component 
in the higher dimensional one is the subset defined by = 0, and consequently also L v = 0, 
which follows from the constraint ( P-16[ ). Again, this implies that the terms involving the link 
7] are not present in the expression ( |3.9| ) for the symplectic potential. Since we also know that 
no other link variables are affected when we insert a new link, we conclude that the symplectic 
structure on the lower dimensional component of the kinematical phase space is equal to that on 
its image in the higher dimensional component. 

The same holds for the mass shell constraints and the Hamiltonian, which is given by ( J3.29| ). If 
any of the transition functions is trivial, thus g x = 1, then it does not contribute to the holonomies 
in ( [3.25D . Once again, it is therefore irrelevant for the definition of the holonomy whether the new 
link has been inserted or not. And the value of the Hamiltonian is consequently also invariant 
under large gauge symmetries. All together, this implies that the quotient space «S/~ is a well 
defined symplectic manifold, and the Hamiltonian, which is actually a linear combination of the 
mass shell constraints, is a well defined function thereon. 

For the free particle system, we had an explicit definition of this manifold. It was the original, 
6n — 4 dimensional phase space spanned by the positions x w and the momenta p w , with the 
restriction to the centre of mass frame imposed. Here we do not have such an explicit definition. 
All we have is a finite atlas, where each chart is a labeled by a graph T, and where the coordinates 
are equivalence classes of kinematical states. Hence, both at the spacetime level and at the phase 
space level we only have local coordinates. A global chart on the phase space can only be 
defined if we also introduce a global coordinate chart on the spacetime manifold [ fRj , |T5| |. But 
then we also have to give up the simple physical interpretation of the phase space coordinates. It 
is not possible to use the relative position vectors and the conjugate transition functions as global 
coordinates on the phase space QT~8|] . 
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4 Einstein gravity 



In this section we are going to derive the phase space structures defined in section [| from the 
Einstein Hilbert action. The starting point is the definition of the particle model at the level of 
general relativity as a field theory on a fixed spacetime manifold Ai. It is convenient to stick 
to the ADM formulation from the very beginning, so that Ai — R X Af is foliated by a space 
manifold Af evolving with respect to some unphysical ADM time coordinate t. Generalizing the 
definition of the previously considered single particle model, we define Af to be a two dimen- 
sional surface, which is either a plane or a sphere, and we cut out with n open discs, so that n 
circular boundaries arise, representing the locations of the particles in space. 

This provides a proper regularization of the conical singularities, and it formally converts the 
matter degrees of freedom associated with the particles into topological degrees of freedom of 
the gravitational field. A set of boundary conditions is imposed to make the boundaries look 
like points in space, and boundary terms are added to the action, to define the coupling of the 
particles to the gravitational field. A comprehensive description can be found in section 2 of [|T3|]. 
For an open universe, we finally impose a fall off condition at spatial infinity. This is a kind of 
asymptotical flatness condition. It defines the centre of mass frame of the universe as a reference 
frame in the sense of [|I9|]. 

To the so defined model we shall then apply the phase space reduction. We set up the Hamilto- 
nian framework, identify the phase space, derive and solve the constraints, and finally the gauge 
degrees of freedom can be divided out. We'll find a finite number of observables, representing 
the physical degrees of freedom of the model, and these are the link variables used in figure |] to 
define the phase space. We shall also recover the symplectic structure and the various kinemat- 
ical constraints defined in the beginning of section 0, which are then shown to be derived from 
the Einstein Hilbert action. 



The space manifold and the particles 

The ADM space manifold Af is a two dimensional, orientable surface with n circular boundaries, 
representing the locations of the particles in space. The global topology of Af is in principle 
arbitrary. But for simplicity, and also for the reasons given in section 0, we shall here restrict to 
two special cases. For an open universe, the space manifold M is a plane M 2 with n open discs 
cut out, as shown in figure 0. For a closed universe, the same n open discs are cut out from a 
sphere S 2 . 

In the first order dreibein formulation of general relativity, the basic field variables are the 
dreibein e M and the spin connection u)^ where /i, u, . . . are formal tangent indices on Ai = 
R x Af. For both open and closed universes, the tangent bundle of the spacetime manifold 
is trivial. Both fields are therefore sl(2)-valued one-forms. We don't have to define them as 
sections in a bundle. The metric and the dreibein determinant are given by 

9fMU = lTi(e^e u ), e = Tr(e^e u e p ) > 0, (4.1) 

where e^ vp is the Levi Civita tensor on Ai. We require the metric to be invertible, and the 
determinant to be positive everywhere in the interior of Ai . What happens at the boundaries will 
be explained in a moment. 
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Figure 9: The space manifold Af for an open universe is a plane R 2 with n open discs cut 
out. For a closed universe, the same open discs are cut out from a sphere S 2 . The circular 
boundary B w represents the location of the particle -k in space. The metric is degenerate 
on the boundaries, so that they look like points in space. The fall off condition at infinity 
requires the existence of a neighbourhood of infinity U C Af, which is embedded into a 
spinning cone. This defines the centre of mass frame and serves as a reference frame. 



In the ADM framework, the fields are decomposed into one-forms and on Af, where 
... are formal tangent indices, and scalars e t and u) t . The induced metric on Af and the 
normal vector are given by 

gtj = ^Triaej), n = \e ij [e h e,]. (4.2) 

For the Levi Civita tensor on Af we have e' lj = e uj . For the foliation to by spacelike, and the 
time orientation of the dreibein to coincide with that of the ADM time coordinate t, the normal 
vector n is required to be negative timelike. A future pointing timelike vector ^ is then mapped 
onto a positive timelike Minkowski vector £ M e M by the dreibein. 

The particle boundaries of Af are denoted by B^. For a circular boundary to look like a point in 
space, its circumference has to vanish. This is the case if and only if the tangent component of the 
dreibein vanishes on the boundary. Let (r, <p) by a polar coordinate system in the neighbourhood 
of the boundary, so that the boundary is at r = 0, and (p increases in counter clockwise direction 
with a period of 2n. The point particle condition then requires 

<Vk = °" ( 4 - 3 ) 

This implies that both the normal vector n and the dreibein determinant e are zero on the bound- 
ary. Therefore, the restriction to negative timelike normal vectors and invertible metrics applies 
in the interior of Af only. For a one-dimensional boundary of the space manifold to look like a 
point, the metric on the boundary has to be degenerate. 

From the spacetime point of view, the two-dimensional cylindrical boundary Rx8 T represents 
a one-dimensional object, the world line of the particle it, and therefore the spacetime metric on 
the boundary has to be degenerate as well. But this is not going to be a problem, because in the 
first order formulation of Einstein gravity in three dimensions, the inverse dreibein never shows 
up. Blowing up the world line to a cylindrical boundary in fact removes all the singularities, 



57 



which are otherwise present in the dreibein and the spin connection describing the gravitational 
field of a point particle [|T3|]. 



Action and boundary terms 

The next step is to define the action functional for the particle model. In the first order dreibein 
formulation, the ADM Lagrangian depends on the fields e M and uj^ and their time derivatives. It 
is the sum of a bulk term, a boundary term for each particle boundary, and for an open universe 
there is also a boundary term at spatial infinity, 

£ = £o + £»-ECA (4-4) 

The bulk term £ is the usual first order Einstein Hilbert Lagrangian, which is given by 

Co = — Jj?x e lj Trfa e j + \u> t T {J + \e t F y ). (4.5) 

The dot denotes the derivative with respect to t, and and are the spatial components of 
the curvature and the torsion, respectively, 

Fij,v = - dyU^ + [w^, u> v ), 

= d^e-v - d v e^ + [u^, e u ] - [u v , e M ]. (4.6) 

This is the same as (3.1) in []T3]], with Newton's constant restored. It is the standard form of the 
Einstein Hilbert Lagrangian in three spacetime dimensions. The variation of the bulk term with 
respect to the dreibein and the spin connection, and with all boundary terms neglected, implies 
that 

= 0, = 0. (4.7) 

These are the vacuum Einstein equations, stating that the spacetime is locally flat outside the 
matter sources. 

The coupling of the particles to the gravitational field is introduced as follows. At each particle 
boundary B w , we add a mass shell constraint C w to the Lagrangian, with a multiplier ( w as a 
coefficient. Formally, the mass shell constraints are given by ( 3.29 ), 

u n - cos(47rGV7v) 
C * ~ 16^ (4 - 8) 

But the holonomy of the particle is now defined as a function of the spin connection on the 
boundary, 

u w = iTr(Pexp f dy^). (4.9) 

The path ordered exponential under the trace represents the Lorentz rotation acting on a vector 
which is transported once around the particle. 

So far, this is a straightforward generalization of the previously defined single particle model. 
The mass shell constraint serves in at this level a two-fold purpose. It depends on the spin 
connection and therefore provides a boundary term for the Lagrangian. One can show that this, 
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together with the point particle condition and the vacuum Einstein equations, implies that the 
particle is moving on a geodesic, and that there is a conical singularity on this geodesic. On the 
other hand, the actually mass shell constraint, which is the equation of motion for the multiplier 
C,r, implies that this geodesic is timelike if the particle is massive and lightlike if it is massless, 
and it fixed the deficit angle of the conical singularity. All this is shown in []T3[], and it now 
applies to each particle independently. 



Asymptotical flatness 

For a closed universe, we now already have a well defined dynamical system. The configuration 
space is spanned by the fields e M and w M , and the multipliers The fields are defined on a 
fixed space manifold J\f, and they are subject to various boundary conditions. We have a well 
defined Lagrangian, which depends on the configuration variables and their time derivatives. The 
boundary term is absent for a closed universe, and the space manifold J\f is compact, so that 



the integral defining the bulk term ( |4.5| ) is finite. 

For an open universe, this is not the case. We have to impose some fall off conditions at infinity, 
and we also have to add an appropriate boundary term C^. Let us require that, far away from 
the particles, the spacetime looks like the gravitational field of a single particle, whose mass M 
is the total energy, and whose spin S is the total angular momentum of the universe. The metric 
at infinity is then that of a spinning cone, which we already encountered in section [ZJ 

ds 2 = -(dT - AGS d0) 2 + dR 2 + (1 - AGM) 2 R 2 d0 2 . (4.10) 

Due to the absence of physical degrees of freedom, or local excitations of the gravitational field, 
we do not have to think about the actual fall off behaviour of the metric at infinity. 

The spinning cone is, up to coordinate transformations, the most general solution to the vac- 
uum Einstein equations in the considered region of spacetime, which admits the definition of a 
centre of mass frame. The condition for a centre of mass frame to exist is the same as for the 
free particle system. The total momentum must be timelike. For the gravitating particles, this 
total momentum is actually the total holonomy. It is defined like the holonomy of the particles, 
but with the spin connection integrated along the circle at infinity. If the total holonomy is not 
timelike, then the universe does not even even admit a proper causal structure [|7|, §]. On the other 
hand, if it is timelike, then it is always possible to choose coordinates so that the metric becomes 

on 



The region where the spinning cone metric applies is a neighbourhood of infinity. A neigh- 
bourhood of infinity in space is a subset hi C A/", which has the topology of K + x S 1 , and whose 
complement is compact. A typical neighbourhood of infinity is shown in figure |^. A neighbour- 
hood of infinity in spacetime is a subset of M. = Ix J\f, so that at each moment of time the 
intersection with the space manifold J\f at that moment of time is a neighbourhood of infinity in 
space. So, the fall off condition to be imposed on the metric is that, at each moment of ADM 
time t, there exists a neighbourhood of infinity hi, and in hi it is possible to introduce a conical 
coordinate system (T, R, <p) so that ( ft.lU[ ) holds. 

Let us translate this into a fall off condition for the dreibein and the spin connection. A possible 
representation of the spinning cone metric is given by 

e M = (d^T + AGS <9 M 0) 70 + <9 M i? 7 (0) + (1 - 4GM) Rd^ {</>), 
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-2GM«V> 70 . 



(4.11) 



One can easily check that this implies = and = 0, so that the spinning cone is indeed 
a solution to the vacuum Einstein equations. However, it is not possible to impose this directly as 
a boundary condition on the basic field variables, because it still involves time derivatives. But it 
turns out to be sufficient to impose only the spatial components of Q4.1 ip as boundary conditions 
at infinity. 

The precise fall off condition is given as follows. At each moment of ADM time t, there 
exists a neighbourhood of infinity hi C A/", a pair of real numbers M and S, and a set of conical 
coordinates (T, R, <p) on hi, so that 

et = (diT + AGSdrf) 7o + <9*#7(0) + (1 - 4GM) Rd^-f'^), 

w< = -2GM d4 lQ . (4.12) 

Note that we do not fix the conical coordinates in any way, or require the conical time T to be 
related in any way to the ADM time T. The conical coordinates are instead considered as addi- 
tional field variables supported on hi, which are related to the dreibein and the spin connection 
by this fall off condition. And it is also useful to note that M and S are implicitly defined as 
functions of the spatial components of the dreibein and the spin connection. 

Now, suppose that the dreibein and the spin connection satisfy this condition at each moment 
of time. Is the bulk term ( ff~5| ) of the Lagrangian then well defined? The last two terms are finite. 
The torsion and the curvature i^j both vanish on hi, and therefore these terms have a compact 
support. What remains is the first term, involving the time derivative. In hi, it is given by 

J- e «Tr(w iCi ) = ¥-e i !d i [(T + AGSi>)d j <j>}. (4.13) 

07TG " Z7T 

Here we added a total time derivative, which simplifies the result slightly. This is of course al- 
lowed, as we may at any time add a total time derivative to the Lagrangian. So, within the neigh- 



bourhood of infinity hi, the integrand in ( |4.5| ) can obviously be written as total spatial derivative. 
To obtain a well defined finite Lagrangian, we just have to subtract the resulting boundary term 
at infinity. Hence, the appropriate term to be added in Q4.4Q is 

£co = -^ I &sd s <p(T + AGS&. (4.14) 

Z7T Joo 

The notation means that the integral is evaluated along a circle at infinity, in counter clockwise 



direction. More precisely, the total Lagrangian is defined as follows. We first integrate ( |43[ ) over 
a connected compact subset A/q C ftf, which is sufficiently large, so that J\f = J\f U hi, and so 
that the boundary cWq is a loop in hi. Then we add the term Q4.14] ), integrated along the boundary 



dj\f . According to ( |4.13[ ), the result is invariant under deformations of the boundary dj\f . Thus 
we may finally also take the limit J\f — ► J\f, where the boundary cWo becomes a circle at infinity. 

With the fall off condition imposed and the boundary term added, the Lagrangian becomes 
finite. However, it is no longer only a function of the dreibein and the spin connection. The 
conical coordinates and their time derivatives show up explicitly in the boundary term. We have 
to include them into the definition of the configuration space. For a given geometry of the ADM 
surface, the conical coordinates are only fixed up to a time translation T i— > T — AT, and a spatial 
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rotation i— *■ + A0. These are the Killing symmetries of the spinning cone. So, we have to 
add two extra degrees of freedom to the configuration space, in order to obtain a well defined 
Lagrangian. 

This is a typical feature of general relativity in the Lagrangian, or Hamiltonian framework 



Q19|]. The additional degrees of freedom represent a reference frame, and the Killing symmetries 
are the possible translations and rotation of the rest of the universe, thus in this case the particles, 
with respect to this reference frame. By definition, this reference frame coincides with the centre 
of mass frame of the universe. It is this feature which is responsible for the various ambiguities 
associated with the reference frame to disappear, when we go over from a single to a multi 
particle model. In []T3[], we found that the boundary term at infinity is not uniquely determined 
by the condition that the Lagrangian should be finite. 

To be precise, this is also not the case here. We are still free to add any function of M and S 
to the Lagrangian. However, one can show that the given boundary term is uniquely fixed by the 
addition condition that the reference frame must not be translated or rotated when time evolves. 
If we derive the equations of motion and take into account that the conical coordinates are now 
also configuration variables, we find in addition to the vacuum Einstein equations the following 
time evolution equations for the conical coordinates, 

e t = (r + 4GS0)7o + i?7(0) + (l-4GM) J R07 / (0), 

u> t = -2GM07O. (4.15) 

Clearly, these are just the time components of 1 1[ ). If we would add a function of M or 
S to the Lagrangian, then there would be additional time translation or rotations of the conical 
coordinates modifying the equations ( |4.15[ ). Physically, this would mean that the reference frame 
is translated or rotated while the universe evolves, and this is of course not what we want. We 
therefore conclude that the Lagrangian is uniquely fixed by the condition that it must be finite, 
and that it has to provide the correct equations of motion for the reference frame. 

Hamiltonian formulation 

Let us now go over from the Lagrangian to the Hamiltonian formulation. This is actually only a 
formal step, because the Lagrangian is already of first order in time derivatives. The phase space 
is spanned by the following field variables on J\f, 

- the spatial components e ; and u?, of the dreibein and the spin connection, 

- the conical coordinates (T, R, <p) defined in some neighbourhood of infinity U. 
They are subject to the following boundary conditions, 

- in the interior of TV, the normal vector n = |e w '[ej, ef\ is negative timelike, 

- on each particle boundary B w , the tangent component of the dreibein vanishes, 

- in the neighbourhood of infinity U, the relation ( |4.12p is satisfied for some M and S. 
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From the total Lagrangian, we can read off the symplectic potential and the Hamiltonian. The 
symplectic potential is the one-form 0, which is given by those terms of C which are linear in 
time derivatives. There is a contribution from the bulk term (fO|), and from the boundary term 



Q4-14D at infinity, 

If Mr 
6 = — — / d 2 xe^Tr(do; i e 7 ) - — / ds d s (f) (dT + AGS d</>). (4.16) 

Here and in the following, the upright letter d always denotes the exterior derivative on the phase 
space, except when it appears immediately behind an integral, where it denotes the measure on 
M or along a curve on M . The Hamiltonian H is defined by those terms on C which do not 
contain time derivatives. There is a contribution from the bulk term, and one contribution from 
each particle boundary, 

U = E ^ - -l^J^x^Tr^tTij + etFij). (4.17) 

This is obviously a linear combination of primary constraints, with the multipliers e t , u> t , and Q 
as coefficients. These are not considered as phase space variables, but rather as free parameters 
entering the Hamiltonian. The kinematical constraints are those associated with the multipliers 
e t and u> t , 

Fij « 0, Tij « 0. (4.18) 
The dynamical constraints are the mass shell constraints, associated with the multipliers Ctt, 

= M ,-cos(4,G m ,) ^ Q 

l07T z G z 

The phase space reduction will now be applied to the kinematical constraints. We are going to 



solve the constraints Q4.18D , and divide out the associated gauge symmetries. The mass shell 
constraints will not be affected by this reduction. It is then already obvious that the reduced 
Hamiltonian is a linear combination of the mass shell constraints, as all other terms in Q4.17P are 



zero when the kinematical constraints are solved. We may therefore in following only focus on 
the symplectic structure. 

Minkowski coordinates 

The general solution to the kinematical constraints ( |4.18| ) is well known. Within any simply 
connected subset of the space manifold, it can be parameterized by a group valued field g : J\f — ► 
SL(2) and a vector field / : M -»• tf 1(2), so that [g, § 

uji = g^dig, e; = g x difg. (4.20) 

It follows that the metric is given by 

9ij = |Tr(e iej ) = iTridJdrf), (4.21) 

and the normal vector becomes 

n = \e ij [ei,ej] = \e ij g'^dj ,d j f]g. (4.22) 
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Figure 10: The triangulated space manifold M is covered by a collection of simply connected 
polygons A € 77. A neighbourhood of infinity hi is divided into simply connected segments 
1/nA, one for each non-compact polygon A € /Too- The links are assumed to be spacetime 
geodesies, and in U the external links 77 € -T^ are spatial half lines in the spinning cone, 
which conical coordinates T — > T n and <fi — > (f> v . 



The field / defines an isometric embedding of the space manifold into Minkowski space. This 
is just a different way to define the Minkowski charts that we already introduced in section 
The components of / are the Minkowski coordinates. The field g represents a Lorentz rotation 
at each point in space. It tells us how the local frame defined by the dreibein at a point in space 
is related to the Minkowski frame defined in the respective chart. 

If we want to parameterize the solutions to the kinematical constraints globally on TV, then 
we have to cover the space manifold by a collection of simply connected subsets. We therefore 
introduce a triangulation, in the very same way as in section [IJ. The only difference is that the 
links A G r are now oriented curves on J\f, and the vertices are the particle boundaries B w . The 
internal links A G r are curves beginning at a particle boundary B 7r _ x , and ending at a particle 
boundary B nx . The external links 77 G extend from a particle boundary B n to infinity. The 
polygons A G II are simply connected, closed subsets of J\f, overlapping along the links. A 
typical such triangulation is shown in figure [TO, 



Now, suppose we are given a field configuration and ojj on Af, satisfying the kinematical 
constraints i^- = and = 0, and a graph r defining a triangulation of J\f. We can then 
introduce a pair of fields g& ■ A — ► SL(2) and /a : A — >• sl(2) in each polygon, so that 

<*>ih = g a dig a, e;| A = g^dif&g^- (4.23) 

The image /a (A) becomes a spacelike surface embedded into Minkowski space. We call this 
the embedded polygon. The embedded polygons are those in figure |3], provided that we impose 
some further boundary condition on the fields g& and /a- 

The first observation is that the fields g& and /a are not uniquely determined by the dreibein 
and the spin connection. The right hand sides of ( |4.23| ) are invariant under the transformations 

g A h-> h A 0a, /a i-> ft-A (/a - «a) ^a\ (4.24) 
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where h A £ SL(2) and «a £ s l(2) are two independent constants for each polygon A E 77. 
This is a Poincare transformation, thus a coordinate transformation in the respective coordinate 
charts. On the embedded polygon /a(A), it acts as a Lorentz rotation by h,A, and a translation 
by t> a- Clearly, these are the redundancy transformations already considered in section 



Link variables 

Referring to the embedding fields g A and /a, we can now also define the link variables. When 
considered as a subset of J\f, a polygon A is not only bounded by the edges A E 7a, but also 



by some segments of the particles boundaries A fl £>?,-. The point particle condition Q4-.3D implies 
that the field /a is constant on these segments, 

^hns* = A'a^JaA'aI^ =0 f A \ B7T = x nA . (4.25) 

The constant x w a is called the position of the particle ix in the polygon A. There are as many 
positions x nA of the particle n as there are polygons A adjacent to n, and these position are 
the corners of the embedded polygons in figure |3|. We also see that the point particle condition 
indeed implies that the particles are pointlike. 

Now, let us assume that all the links are geodesies in spacetime. This can be regarded as a 
gauge condition in general relativity, restricting the possible foliations and also the embeddings 
of the triangulated space manifold J\f into the spacetime. If this is the case, then the image /a (i]) 
of each link is a straight line in Minkowski space. Let us introduce on A an affine parameter s. 
For internal links A E To, it is chosen so that s E [0, 1], and A(0) = 7t_a and A(l) = ti\. For 
external links, we choose s E [0,infty), and for reversed external links s E (— oo,0], and in 
these cases s measures the proper length. Then we define 

z r , = d s f Av , Tj E r . (4.26) 

The vector z\ is a constant, because the link A is a geodesic. For external links, it is a spacelike 
unit vector. For internal links, it is the relative position vector, 

z\ = x nxAx - x n _ xAx , A E r . (4.27) 

The transition functions are defined as follows. Consider a pair of links ±77, and the adjacent 
polygons A± n . The both links represent the same curve on TV, and therefore the dreibein and the 
spin connection are equal on A and —A. But the fields g Ax and /a a on A, and <7a_ a and /a_ a 
on —A are in general different. On the other hand, we have seen above that if two embeddings 
define the same dreibein and the same spin connection, then they are related by a Poincare trans- 
formation of the form ( [4.24] ). Hence, for each pair of links ±77 there exists a constant g\ E SL(2) 
and a constant fx E sl(2) so that 

c/a_ a = g\gA x , fA x = gx (/a a - fx) g\\ Xe r. (4.28) 

Exchanging A A and A_ A in the first equation, and differentiating the second equation with re- 
spect to the curve parameter s, where we have to take into account that the reversed link is 
oriented into the opposite direction, we derive the following relations, 

g x = gx 1 , z-x = -gx z x gx l , \e r. (4.29) 
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Hence g\ is the transition function, and the relation ( |3.3D is recovered for both internal and 
external links. The translational component fx of the transition function has never shown up in 
section because there are link variables that refer to the absolute positions of the polygon in 
the embedding Minkowski space. 

The conical frame 

There is one more gauge condition that we have to impose. The external links must be spatial 
half lines in the spinning cone. Hence, if we follow a curve i] G on J\f to infinity, then at some 
point the curves enters the neighbourhood of infinity U, and in the limit the conical coordinates 
are given by 

T — > T v , 4> -> <f> v , R -> oo, r\ G r^. (4.30) 

Moreover, there is also a further condition to be imposed on the embedding of the non-compact 
polygons. Let us choose the neighbourhood of infinity U so that, for each non-compact polygon 
A G i/oo, the intersection U n A is a simply connected segment, as indicated in figure [TTJ This is 
always possible, because given any neighbourhood of infinity we can always go over to a smaller 
one which has this property. 

The segment U PI A is then not only embedded into the spinning cone, but also into Minkowski 
space. The first embedding is provided by the conical coordinates (T, R, <p), the second by the 
field /a- And we still have the freedom to apply a Poincare transformation to f A . Using this, 
we can always achieve that the two embeddings are related by a local isometry of the form 
( J2. 19| ). There is then the following relation between the conical coordinates (T, R, 0), and the 



Minkowski coordinates f A in U n A, 

/a = (T + r A (0)) 70 + #7(0 - «a(0)), (4.31) 
where ta and a A are two linear functions satisfying 

r A (0) = AGS, a' A (<t>) = AGM. (4.32) 
For the fall off condition ( [4. 12[ ) to be satisfied, the field g A is then necessarily given by 

9A = e -«AWW2. ( 43 3) 

Using this, we can finally also derive the special transition functions g± v , and the unit vectors 
z± v for the external links rj G r^. For this purpose, we have to evaluate the fields g A and f A on 
the external links ±77, in the limit R — ► 00. What we find is 

SA_„ - e^ 2 , f A _ v - (T„ + rA_„(0„)) 7o + #7(0, + 0"), 
SA„ - e^°/ 2 , /a, - (T„ + r A ,(0,)) 70 + #7(0, + ^ + ), (4-34) 



where we introduced the deviations already defined in ( |2.24| ), thus 



9+ = -aA v (<f> v ), V = -«a_„(0,)- ( 4 -35) 



65 



Now, consider first the vectors z± v , as defined in fl4.27| ). In the limit R — > 0, the conical coordi- 
nate R asymptotically defines the proper length on a spatial half line, and therefore we can use 
s = ±R as a parameter on the links ±r/. It then follows from ( |4.34{ ) that 



7(0^ + e+ 



-7( 



e: 



v e r 



(4.36) 



which is the same as (|3.5[). To find the transition functions g± v , we have to use the defining 
relation ( |4.28| ), which can also be evaluated in the limit R — ► 0. What we get is 



9n 



-1 

9a_„ gA v 



-{e+-d-)i /2 



We can therefore make the definition 



87rGM„ = 9t - 9; 



9n 



-4nGM v -yo 



(4.37) 



(4.38) 



and recover the definition ( |3T4"| ) of the energy M v , and also the second equation in (p.6|). The first 
equation in ( |3.6| ) follows directly from the definition fl4.35| ). With A = A^ = A_^/ we have 



9 V , - 8+ 



-a A [ 



+ a A ( 



-4GM* 



(4.39) 



Finally, what remains is the average condition fl3.7|). In section [Z] we have seen that this is equiv- 
alent to the condition that the integral of the functions a(4>) over the circle at infinity vanishes, 
hence 

/ dsd s (f)a((f)) =Y1 I ds9 s 0a A (0) = 0. (4.40) 

ioo a , — rT JoonA 



This is the same as ( |2.28| ), we only changed the notation slightly, because now we have consider 
this as an integral along a curve on M, and is a field variable onW C ftf. The notation oo fl A 
means that the integral is evaluated along the segment of the circle at infinity which belongs to 
the non-compact polygon A. 



The symplectic structure 

With these formulas at hand, we can now compute the symplectic structure on the subspace 



defined by the kinematical constraints fl4.18p . The embedding fields g A and /a can be used as 



coordinates on this subspace. To derive the symplectic potential, we insert the parameterization 
( ]4.23[ ) into the original expression ( )4.16D , 



1 r Mr 
= —— d 2 x^Tr(du> iej ) - — / ds0 s 0(dT + 4GSd0). (4.41) 

07TG JM ' lH Joo 

For simplicity, let us first ignore all boundary terms at infinity, or consider a closed universe 
where these are absent. The symplectic potential is then given by the bulk term only. To write 
this as a function of the embedding fields g A and f A , we first have to split the integral into a sum 
of integrals over the individual polygons, 

o = 77^7 E /d 2 x^Tr(d^e,) = -i-E [ d 2 x d t Tr(dg A gK 1 d 3 f A ). (4.42) 

07T(jr » n JA 07T(jr . „ JA 
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Here we have inserted < J4.23| ), and performed some simple algebraic manipulations. The inte- 
grand is thus a total derivative, and we can write the result as a boundary integral, 



o = E / dsTi(dg A g A 1 d s f A ). 
o7rG £r n J OA 



(4.43) 



The boundary <9A of the polygon A is always traversed in counter clockwise direction. It consists 
of the edges A G r A , and the segments of the particles boundaries belonging to the polygon A. 
The particle boundaries, however, do not contribute to the integral, because according to the point 
particle condition f A is constant, and therefore the integrand vanishes. What remains is 



o = 77^E E I dsTr(dg A g A 1 d s f A ). 

07TG TTn \Tt 



(4.44) 



Now, we may equally well sum over all links A e F, because each link appears exactly once as 
a boundary of a polygon, 



= E J x dsTr(dc/ AA g Ax dJ Ax ). 



(4.45) 



For a closed universe, all links are internal links, thus r = T . Using this and the decomposition 
r = r + n JL, we can also write 

I p 

AG-T+ 

We can then use the defining relations Q4.28D of the transition functions, to express g A _ x and 
/ A _ A in terms of g Ax and f A . This gives 



(4.47) 



AG/ + 



Now, the transition function g\ is a constant, and therefore the integrand is again a total deriva- 
tive. Using this and the definition ( |4.27| ) of the relative position vector Z\, we get 

1 



6>n 



8ttG 



E Tr(gi A dg x z x ) 



(4.48) 



Aer+ 



So, for a closed universe we have shown that the symplectic potential, derived from the Einstein 
Hilbert action, is indeed given by the expression ([3 .9Q. The terms associated with the external 
links are in this case not present. 

For an open universe, we perform the same calculation, but we also have to take into account 
some boundary terms at infinity, and the special properties of the external links. Let us first 
consider the external links, still ignoring the boundary terms at infinity. The total symplectic 
potential can then still be written as ( |4.45| ), but in ( |4.46| ) and consequently also in ( 4-.47D we have 
to sum over all external links as well. We shall now show, however, that the external links do not 
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contribute to this sum. Consider an external link rj e i^, and the corresponding term in ( |4.47| ), 
which is given by 

-J- [ds Tr(g^ 1 dg v d s f Av ). (4.49) 
The integrand is still a total derivative, but the integral cannot be expressed like (|4.48|), because 



a relative position vector assigned to the external link r\ does not exist. However, the integrand is 
actually zero. This can be seen as follows. According to ( |4.37[ ), the transition function g n is an 



element of the subgroup S0(2), which implies that the term g^ 1 dg v is proportional to 70. On 
the other hand, according to ( |4.26[ ) and ( |4-.36| ), the derivative d s f Ari is orthogonal to 70, and thus 
the integrand is zero. 

So, there is no contribution to the symplectic potential from the external links. The only 
additional terms for an open universe are the boundary terms at spatial infinity. There are two 
types of such boundary terms. Let us first consider the boundary term in ( |4-.41| ). It is an integral 
along the circle at infinity. We split it into segments belonging to the external polygons A £ 11^, 

--— V / dsAGM8 s <P{dT + AGSd<P). (4.50) 
8ttG Ae77oo JoonA 

Remember that the abbreviated notation ooflA means that the integral in evaluated along the 
circle segment which defines the boundary at infinity of the non-compact polygon A e iloo- The 



integrand is written in this particular form, because according to ( |4.32[ ) we have 

AGMd s (j) = d s a A ((j)). (4.51) 

Hence, ( |4.50| ) is equal to 



~TTf £ /d S 9 s a A (0)(dT + 4G5d0). (4.52) 
87rG Aen x JoonA 

There is then a second boundary term at spatial infinity, which arises from the boundary integral 



Q4.43| ). The boundary of each non-compact polygon also includes a segment of the circle at infin- 
ity, which we did not take into account so far. Hence, yet another contribution to the symplectic 
potential is 

l — I dsTi(dg A g A 1 d s f A ). (4.53) 

57rG AeI7 a JoonA 



The path along which this integral is evaluated lies entirely within the segment U PI A of the 
neighbourhood of infinity. We can therefore insert the expressions ( |4.3 1| ) for f A and Q4-.33D for 



g A . Using this, one finds that Q4.53D is equal to 

1 

8^G 



X — /dsd(a(0)) {d s T + AGSd s cj>). (4.54) 

Let us add to this a total exterior derivative, which we are always allowed to, as the only relevant 
object is the symplectic two-form i? = dO. Then we have 

V / dsa A (4>)[d s (dT + 4GSd<f))+4GdSd s 4>}. (4.55) 
8vrG Aei7oo JoonA 
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The term proportional to dS is the expression which is set to zero in Q4-.40| ). It therefore vanishes 



What remains is, together with fl4.52| ), again a total derivative, 



A€77 t 

The integral can be evaluated and becomes 
1 



0oo = -^ E / dsd 8 [a A (<l>)(dT + 4GSd(l>)]. (4.56) 

07TOr JoonA 



87rG Aei7 c 



E K(<M (dT v , + AGStyrf) - «a(^) (dT„ + 4GS , d0„)]. (4.57) 



Here we used that the segment of the circle at infinity belonging to the polygon A is bounded 
by the two external link r\ E and — r] E -T-oo* where A = = A.- v '. And we inserted the 
conical coordinates on these edges, which are given by flOOp . Finally, we can use the definitions 



Q4.35D and rearrange the summation, which gives 



Ooo = E (tf - d v) ( dT ^ + 4GSd0„). (4.58) 



And now, we just have to use ( |4.38| ) and make another definition, namely 

L V = -AGSM V , (4.59) 

to obtain the final result. The contribution to the symplectic potential from the boundary at 
infinity is, again up to a total derivative, 

Ooo = -J2 M v (dT v + AGS d<j> v ) = E ( T v dM v + L v d<f) v ). (4.60) 



This completes the technical part of the derivation. The sum = 6^ + Oq is the same as fl3.9j). 
Now we just have to sort out what this equality actually means, and how it fits together with the 
definitions made in section |3|. First of all, it follows that all those variations of the embedding 
fields <7a and /a are gauge symmetries of general relativity, which preserve the values of link 
variables showing up in ( |4.48p and fl4.60p . In other words, the link variables are the observables 



of the particle model, in the sense that all physical degrees of freedom are encoded in these 
variables. This is what we had to assume in the beginning of section |3|. 

Let us consider the quotient space of all field configurations g& and /a, divided by the gauge 
transformations that preserve the values of the link variables. This quotient space is obviously 
the kinematical phase space Sp defined in section ^ To see this, we first note that by definition 
the observables, hence the link variables are the coordinates on this space. And furthermore, 
they are not independent but subject to the kinematical constraints ( |3.14| - pTT6| ). It is therefore not 
the extended phase space Q r which is recovered here, but the kinematical subspace Sp C Qr- 
To check this, we verify that the constraints defined in section |3| are indeed satisfied by the link 
variables introduced above. 

The constraint ~ for each compact polygon A E n is an immediate consequence of 



the definition ( |4.27D . The same applies to the constraint for each non-compact polygon 



A E iloo. This follows from the definition of the clocks in ( |O0[ ), the fall off condition ( |4.31| ) 



for /a at infinity, and again the definition ( |4.27| ) of the vectors Z\. And finally, the constraint 



J v w for each external link rj E simply follows from the definition ( [4-.59P of the auxiliary 
variable L v . So, everything fits together with the definitions made in section |3|. 
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Appendix 

Here we summarize some facts about the spinor representation of the three dimensional Lorentz 
algebra sl(2) of traceless 2x2 matrices, and the associated Lie group SL(2). A more compre- 
hensive overview, using the same notion, is given in Jl3p. As a vector space, sl(2) is isometric 
to three dimensional Minkowski space. An orthonormal basis is given by the gamma matrices 

7o=(_° 1 J), 7i= (J I), 12= (J (A.1) 

They satisfy 

lalb = Vab 1 - Zabc 7°) ( A - 2 ) 

where a,b = 0, 1, 2. The metric r] ab has signature (—,+,+), and for the Levi Civita symbol 
£ abc we j^yg £ o\2 _ -L Expanding a matrix in terms of the gamma matrices, we obtain an 
isomorphism of si (2) and Minkowski space, 

v = v a la v a = |Tr(«7 a ). (A.3) 

Some useful relations are that the scalar product of two vectors is equal to the trace norm of the 
corresponding matrices, and the vector product is essentially given by the matrix commutator, 

|Tr (vw) = v a w a , \ [«, w] = -e ahc v a w b 7c . (A.4) 

Sometimes it is useful to introduce cylindrical coordinates in Minkowski space, writing 

v = Tj + pj(ip), (A.5) 

where r is the time, p > the radial, and (p the angular coordinate, which is redundant for p = 
and has a period of 2n. The unit vector y(ip) defines the angular direction if in Minkowski space. 
This vector and its derivative "/((p) are two orthogonal unit spacelike vectors, forming together 
with 70 a rotated orthonormal basis, 

~Y(<p) = cosy? 71 + sin tp *y 2 , j'(^p) — cos ip~f 2 — smip-ji. (A. 6) 
Useful relations are 

7o7(^) = l'(<f), 7o7'(^) = -70), (A.7) 

and 

7(<pi) 7(^2) = cos(v9i -ip 2 )l + sin(v9! - <p 2 ) j - (A.8) 
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The Lie group SL(2) consists of matrices u with unit determinant. The group acts on the algebra 
in the adjoint representation, so that 

v i— > u l vu. (A. 9) 

This defines a proper Lorentz rotation of the vector v. The conjugacy classes of the algebra 
are characterized by the invariant length v a v a = |Tr(-u 2 ). For timelike and lightlike vectors 
v (with |Tr(t; 2 ) < 0) we distinguish between positive (v° = |Tr(«T ) > 0) and negative 
(v° = |Tr(t;7 ) < 0) timelike and lightlike vectors, respectively. Special group elements are 
those representing spatial rotations about the 7 -axis. The matrix 

u = e^ o/2 = cos(cj/2) 1 - sin(cj/2) To, (A. 10) 

represents a rotation by uo in counter clockwise direction, 

7o h-> w _1 To u = To, tM h-> u _1 t(</?) u = t(v? + u). (A.ll) 

A generic element u e SL(2) can be expanded in terms of the unit and the gamma matrices, 
defining a scalar u and a vector p e sl(2), 

U = Ml+ P a Ta =>• V=V a ia- (A.12) 

The scalar m is half of the trace of n, and the vector p is called the projection of The determi- 
nant condition implies that 

u 2 =p a Pa + l. (A.13) 

According to the property of the vector p, we distinguish between timelike, lightlike and space- 
like group elements u. For a timelike element we have — 1 < u < 1, and it represents a rotation 
about some timelike axis, which is specified by the vector p. The angle of rotation uj obeys 
u = cos(c<j/2). Lightlike elements are null rotations, and spacelike elements are boosts. 

Finally, a parity transformation, which acts on Minkowski space as a reflection of space, can 
be defined as follows, 

v h-> -T2 vj 2 - (A. 14) 

It represents a reflection of the T2-axis, 

To ^ To, Ti | - > 7i> 72 !-> -T2- (A.15) 
Note that this is not a proper Lorentz rotation, because T2 is not an element of the group SL(2). 
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